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Chapter  I.  INTRODUCTION 


1.  Outline 


We  shall  presen^^unified  treatment  of  boundary  vaiue  prob¬ 
lems  involving  a  system  of  ionized  particles  (i.e.  a  plasma)  whose 
singlet  phase- space  density  is  presumed  to  satisfy  the  linearized 

IS  ^ 


form  of  the  Landau-Vlasov  equation^  Sperifira.lly, -we  shall  make 

4  normal  mode  expansion  in  the  singular  eigenfunctions  :>f  the 
coupled  Maxwell  and  Landau-Vlasov  equations^  ¥  ■-  j  >  , 

This  first  chapter,  by  way  of  introduction,  contains  a  descripX 
tion  of  the  nonstatistical  approach  to  plasma  oscillations,  followed  by 
an  outline  of  the  ideas  of  various  writers  who  have  derived  the  Vlasov 
equation  using  a  variety  of  approaches.  In  the  last  section  of  this 
chapter  we  explicitly  present  the  basic  equations  to  be  considered  in 
the  remainder  of  the  text. 

In  chapter  II  we  exhibit  normal  mode  solutions  of  the  coupled 
equations  for  the  case  of  transverse  waves  of  fixed  frequency.  The 
spectrum  and  orthogonality  properties  of  these  modes  are  discussed 
and  a  fundamental  completeness  theorem  is  proved. 

In  chapter  III  we  treat  the  problem  of  reflection  of  electromag¬ 
netic  radiation  from  a  plasma  half  space  and  slab,  and  discuss  the 
modifications  necessary  when  the  plasma  is  located  in  a  steady-state 
magnetic  field. 
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Chapter  IV  includes  an  assortment  of  topics.  The  initial 
value  problem  for  transverse  modes  is  considered  from  the  stand¬ 
point  of  both  the  singular  eigenfunction  expansion  and  the  more 
conventional  Laplace  transform  treatment.  Since  this  is  the  most 
familiar  problem  is  plasma  oscillations,  only  the  results  are  pre¬ 
sented  and  these  only  in  outline  form.  The  long-time  "Landau"  damping 
that  results  is  discussed  along  with  the  corresponding  damping  in  the 
transverse  mode  initial  value  problem.  The  normal  modes  for  fixed 
frequency  longitudinal  oscillations  are  derived  and  the  corresponding¬ 
completeness  theorem  is  proved.  The  eigenfunction  expansion  is  then 
used  to  solve  longitudinal  mode  boundary  value  problems,  resulting  in 
expressions  for  the  penetration  of  an  oscillatory  longitudinal  electric 
field  into  a  plasma  half-space  and  the  impedance  of  a  plasma -filled 
parallel  plate  capacitor.  Finally,  we  mention  a  rather  indirect  applica¬ 
tion  of  the  singular  eigenfunction  method  to  the  problem  of  electron 
migration  in  a  discharge  tube. 

2.  The  Plasma  as  Dielectric 

The  most  well  known  theory  of  plasma  oscillations  is  Langmuir's 
characterization  of  a  plasma  as  a  collection  of  mobile  electrons  at  zero 
temperature,  and  essentially  immobile  ions.*  According  to  this  very 


In  fact,  Langmuir*  was  the  first  to  use  the  term  "plasma"  to  denote 
a  region  in  an  ionized  gas  containing  equal  numbers  of  ions  and  electrons. 
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rudimentary  model,  an  electron,  normally  at  rest,  moves  under  the 
influence  of  the  electric  field  which  arises  from  the  displacement  of 
all  other  electrons  from  their  equilibrium  positions.  From  this 
interdependence  of  electric  field  and  electron  positions  originates  the 
characteristic  collective  behavior  of  the  plasma.  In  particular,  in 
this  simple  model  the  plasma  behaves  like  a  homogeneous  dielectric 
medium  (with  respect  to  sinusoidal  disturbances  of  frequency  0)  ) 
with  permeability 


(1.1) 


e  and  m  are  the  electron's  charge  and  mass  respectively,  and  n 
is  the  electron  density. 

It  follows  that  the  plasma  can  support  longitudinal  oscillations •• 
"plasma  oscillations  "--of  the  frequency 


(1.2) 


for  which  the  permeability  vanishes.  That  is,  Maxwell's  equations 


r*  -a 

VxE 


7>t 


lH 

c  T,t: 


£  V*  E  «  V' B  *  0 


(1.3) 


are  satisfied  for  disturbances  of  the  form 
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(1.4) 


for  any  wave  vector  k  . 

These  oscillations  can  never  be  utilized  to  propagate  a  signal 
(the  group  velocity  9(J^/9k  vanishes  identically),  but  a  plasma  can 
be  caused  to  oscillate  in  such  a  manner  that  surfaces  of  constant 
phase  will  move  forward  -  "a  situation  like  the  familiar  barber  poles 
which  appear  to  move  steadily  without  rising"^ 

The  characterization  of  a  plasma  as  a  dielectric  medium  whose 
dielectric  properties  are  determined  by  relating  electron  positions  to 
the  electromagnetic  field  and  vice-versa,  places  plasma  problems  on 
an  equal  footing  with  problems  of  the  electromagnetic  behavior  of  more 
familiar  materials.  In  this  approach,  the  plasma  is  characterized 
entirely  by  its  permeability,  and  its  unique  behavior  stems  from  the 
peculiar  frequency  dependence  of  the  permeability. 

For  example,  consider  the  case  of  an  electromagnetic  wave 
normally  incident  on  a  plasma  half  space.  The  familiar  formulas  of 
electromagnetic  theory  yield 


(the  index  of  refraction)  (1.  5) 


and 


T1  (m)  =  f£(cO) 


-5 


T(“) 


•n(<o) 

-  (the  transmission  coefficient)  (1.6) 

1 1  +n((o)l‘ 


0  u; «  tOp 

A  complete  exposition  of  the  applications  of  this  approach  to 

problems  of  wave  propagation  in  ionized  atmospheres  may  be  found 

3  4 

in  the  review  article  by  Mimno  .  Also,  Ford  exhibits  the  various 
modes  of  oscillation  of  a  zero  temperature  plasma  under  the  influence 
of  a  constant  magnetic  field,  and  utilizes  them  in  the  solution  of 
boundary  value  problems  in  the  transmission  and  reflection  of  electro* 
magnetic  radiation  from  a  plasma  half  space. 


3,  The  Vlasov  Equation 

A  more  nearly  rigorous  description  of  the  plasma  as  a  collec¬ 
tion  of  charged  particles  with  random  motions  demands  a  statistical 
treatment.  This  is  given  through  the  Vlasov,  or  "collisionless- 
Boltzmann"  equation. 

Let  (r,  u,  t)  denote  the  phase  space  single  particle  density. 

5 

It  obeys  the  Boltzmann  equation 


II 

dt 


eo/ftsiw 


(1.7) 
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where  a  is  the  local  acceleration  produced  by  external  forces 
(if  any)  and  the  right-hand  side  may  be  evaluated  explicitly  by 
examining  the  binary  collision  mechanisms  in  detail.  Vlasov^ 
noted  that  if  the  plasma  is  sufficiently  ionized  and  not  too  dense  (as 
is  the  case,  for  example,  in  the  ionosphere),  the  individual  electron 
collisions  via  both  short  range  and  Coulomb  interactions  occur  with 
frequencies  far  lower  than  the  plasma  frequency.  On  the  other  hand, 
the  long-range  effects  of  the  electromagnetic  forces  due  to  every  other 
charged  particle  ("many-particle"  collisions)  nuiy  be  taken  into  account 
through  the  a  *  V  term.  Thus  Vlasov  asserted: 


'?!)  s  0- 


a  = 


(1.8) 


(1.9) 


where  the  plasma  particles  themselves  serve  as  the  source  of  the 
electromagnetic  fields  E  and  B  .  Hence  the  electromagnetic  fields 
are  made  self-consistent  through  the  coupling  of  the  Vlasov  equation 


(1.10) 


and  Maxwell's  equations. 


4.  The  Vlasov  Equation  (continued) 

Vlasov's  heuristic  "derivation",  despite  its  aesthetic  appeal, 
is  actually  quite  unrealistic.  The  Boltzmann  equation  rests  on  the 
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hypothesis  that  the  particles  spend  a  negligible  fraction  of  their  time 
colliding,  whereas  the  electrons  are  constantly  interacting  through 
the  Coulomb  force.  On  the  other  hand,  the  separation  of  Coulomb 
effects  into  binary  collisions  and  an  "external"  interaction  via  the  sel  - 
consistent  field  possesses  no  rigorous  meaning.  In  actual  fact,  the 
"smearing"  of  an  electric  field  produced  by  a  collection  of  point 
charges  into  a  smoothly  varying  self-consistent  field  must  be  viewed 
as  some  sort  of  an  approximation. 

The  nature  of  this  approximation  appears  directly  from  con¬ 
sideration  of  the  Liouville  equation  for  an  N  particle  system  with 
binary  interactions: 


N 


(1.11) 


where  </)  is  the  binary  interaction  potential. 

Integration  of  this  master  equation  over  position  and  momentum 
coordinates  of  all  but  s  particles  {s=l,  2,  ...  N)  leads  to  a 
hierarchy  of  equations,  the  familiar  B-B-G-K-Y  hierarchy.  If  we 
define  reduced  distribution  functions 


. -O.V)  =vMT„t 

isSti 


(1.12) 


e 


X 


and  set 


4(li- -■?)!)  = 


(1.13) 
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we  obtain  the  equation  for  • 

and  similar  equations  expressing  in  terms  of  7g^.i  ■ 


14) 


Ti%  A  ^ 

Ts 


T  f  = 

^  “.i 


iS| 


.  e 

m  ^  ''i  lr-f,l 


(1.15) 


In  order  to  close  the  hierarchy,  some  statistical  assumption 
is  necessary.  In  particular,  if  we  assume  that  there  are  no  correla¬ 
tions  between  particle  motions,  then  the  two -particle  probability 
distribution  is  the  product  of  the  respective  single -particle  probability 
distributions 


(1.16) 


and  the 


equation  decouples  from  the  rest  of  the  hierarchy; 


l-r-Ji-l 


(1.17) 


This  is  just  Vlasov's  equation  with  the  self-consistent  field  obtained 

7 

through  the  solution  of  Poisson's  equation.  Harris  ,  starting  from 
Liouville's  equation  with  the  correct  (retarded)  electromagnetic  inter- 
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action,  instead  of  (1.13),  shows  that  the  resulting  Vlasov  equation 
contains  the  Lorentz  force  term  as  well,  as  might  be  expected. 

Intuitively,  the  assumption  of  no  correlations  breaks  down 
is  two  particles  spend  much  time  "close"  to  each  other.  Hence, 
the  factorization  (1.16)  rests  on  the  assumption  that  the  density  is 
low  enough  so  that  individual  particle  effects  may  be  neglected  with 
respect  to  long-range  effects.  The  criterion  for  this  may  be  taken 
to  be 

nl*»l  <‘>8' 


where 


\ 


'k 


the  Debye  shielding  length.  (1.19) 


Equation  (1.18)  implies  that  the  long-range  effects  predomiiuite,  since 

a)  it  insures  that  many  particles  lie  within  the  range  of  the 
shielded  potential  of  a  single  electron,  and 

b)  by  virtue  of  the  relation 


it  has  the  effect  that  the  average  potential  energy  of  a  binary 
interaction  is  much  less  than  the  average  kinetic  energy  of 
a  particle. 

Laboratory  plasmas  and  plasmas  of  astrophysical  interest  (solar  corona, 


gaseous  nebulae,  etc.)  satisfy  the  criterion  (1.18)  with  several  orders  of 
magnitude  to  spare. 


10 


g 

Rosenbluth  and  Rostocker  have  developed  a  consistent 

3 

expansion  of  the  hierarchy  (LIS)  in  powers  of  g  =  l/n 
similar  to  the  Mayer  cluster  expansion.  Write: 

It  +  Z  m^i) 

•  ^  •  s.  ' 


^ZT(j,k,l)|r'3JU)  +ZP0WP(Wf 


(1.21) 


*1 


where  the  pair,  triplet . correlations  are  defined  recursively. 

For  example: 

Then  if  dimensionless  units  are  chosen  in  terms  of  the  basic  length 

Xp.  and  the  basic  time  1/ U  ,  (1.15)  becomes: 

u  p 


3  +  7-u  4r  -  2.7 


IS/ 


=  \ka\. 

C«f  ''  ^ 


(1.23) 


and  one  may  assume  an  expansion  of  the  form 


T,  .^‘•^'5'^^'*'+  ... 


(1.24) 
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where  'Jr  order  g"  .  To  order  g®,  the  third  term 


in 


(1.23)  may  be  neglected  and  a  solution  is.  seen  to  be: 

5 


7^ = I  (1.^5, 

where  7^°^  satisfies  the  Vlasov  equation. 

Although  this  is  enough  for  our  purposes,  it  should  be  noted 
that  the  expansion  may  be  carried  to  higher  orders  in  g  quite 
straightforwardly.  In  particular,  Rosenbluth  and  Rostocker  obtain, 

but  do  not  solve,  the  integro-differential  equation  for  P  in  terms  of 

m  1 

,  correct  to  order  g  . 

9 

Bogoliubov  treats  the  deeper  question  of  whether  or  not  a 
kinetic  description  is  even  possible.  That  is,  under  what  circumstances 
does  the  single  particle  distribution  obey  an  equation  of  the  form 


(1.26) 


In  this  notation,  A(r,  u  |  T^)  is  an  expression  which  depends  upon 
time  functionally  through  8  is  completely  deter¬ 

mined  by  the  distribution  *7*^  at  time  t. 

Bogoliubov  notes  that  after  a  time  t^  after  the  system  is  set 
into  motion  the  correlations  become  synchronized  with  the  single 
particle  distributions  to  the  extent  that  the  s-particle  distributions 
depend  upon  time  only  through 
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7,(t, . -a,  I?:) 

Of  course,  for  long  times,  the  particle  motions  become  entirely 
uncorrelated. 


%  T - ►  I 

^  t— •»  ic,  ' 


(1.28) 


an  expression  which  is  indeed  of  the  form  (1.27).  However,  Bogoliubov 
argues  that  this  process  of  synchronization  takes  place  very  rapidly, 
in  the  time  it  takes  an  electron  to  travel  a  Debye  length 


_L 


to. 


and  it  is  from  this  time  onward  that  a  kinetic  description  exists. 
Granting  this  assumption,  one  then  expands 


(1.29) 


il5,0  =  . + 


(1.  30) 


^3 

7)t 


(1.  31) 


and  requires  that  (1.23)  be  satisfied  to  all  orders  in  g.  This  results 
in  a  set  of  differential  equations  for  the  functionals  subject  to 

the  boundary  condition  (1.28).  The  zeroth  order  result  is  again 


(1.32) 


along  with  the  Vlasov  equation. 
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Guernsey  extends  Bogoliubov's  treatment  to  a  detailed  study 
of  the  approximation  to  order  .  He  obtains  an  integral  equation 
for  7 terms  of  ,  and  solves  it  by  means  of  the  technique 
described  in  Chapter  II,  section  5  of  this  work.  The  result  when  sub¬ 
stituted  into  (1.14)  yields  a  correction  to  the  Vlasov  equation  in  the 
form  of  an  additional  collision-like  term.  In  fact,  he  shows  that  the 
presence  of  this  collision  term  insures  that  an  H-theorem  is  satisfied 
and  the  system  is  driven  to  equilibrium.  I.e. 


i  ?  0  ,1.33, 

with  equality  only  when  is  Maxwellian. 

11  12 

Balescu  and  Prigogine  ’  obtain  the  same  results  as  Guernsey 
by  means  of  a  direct  expansion  of  the  full  Liouville  equation  using  dia¬ 
gram  techniques  to  evaluate  the  complicated  sums  that  arise. 


5.  Linearization 

Our  starting  point  then  is  the  Vlasov  equation  which  represents 
the  effect  of  the  Lorentz  force  on  the  electron  distribution  function, 


and  Maxwell's  equations,  with  the  charges  and  currents  being  given 
explicitly  in  terms  of  the  velocity  moments  of  the  electron  distribution. 
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VxE 


c  -jr 


i4rce 


-N. 


-A  ^ 

Vx5 


7'B  -  0 


(1.35) 


We  have  assumed  that  the  ions,  by  virtue  of  their  greater  mass,  are 
relatively  immobile  and  simply  comprise  a  uniform  background  of 
density  .  This  restriction  is  relaxed  in  chapter  IV.  These  equa¬ 
tions  may  be  linearized  under  the  assumption  that  the  electron  distri¬ 
bution  does  not  deviate  substantially  from  its  spatially  uniform  steady- 
state  value,  l.e. 

(  n  is  the  electron  density) 

f  Cr,  "u,  t)  is  to  be  regarded  as  a  "small"  quantity. 

Since  n  =  for  a  neutral  plasma,  and  since  we  will  be  deal¬ 
ing  with  situations  in  which  the  steady-state  electron  distribution  does 
not  contain  currents,  the  charge -current  distribution  will  be  due  entirely 
to  f  (r,  u,  t)  . 

(1-37) 


(1.38) 
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Thus  the  electric  and  magnetic  fields  are  also  "small"  quantities 
and  to  first  order  in  small  quantities  we  obtain. 


“  m  (1.39a) 


-A  ^ 

V-E 


1“^ 
c  ‘if 


=  4TCe|fcl^u 


V’-6  -0 


(1.  39b) 


(1.39c) 


6,  Separation  into  Longitudinal  and  Transverse  Modes 
It  is  easy  to  see  that  for  isotropic  9^  (u  )  ,  i.e. 

7,(u)  =  (1.40) 

(the  Maxwell-Boltzmann  distribution  is  such  a  function),  the  equations 
(1.  39)  may  be  separated  into  three  independent  sets  of  equations.  First 
note  that  in  view  of  (1. 40) 

(uxi)-K7,  =  T.aSDj  =  0  ,1.«) 

We  shall  seek  plane  wave  solutions  of  the  coupled  equations  (1.  39). 

Set  d/dx  =  d/dy  -  0  ,  write  u  =  (u  ,  u  ,  u)  and  define 

X  y 


(1.42) 


(1.43) 


Now  for  isotropic 


7 


o 


(no  sum  convention  implied) 


where 

0^)  *  i”  !  "  -J"!  =  -  FW  "• «' 

''-0#  ''-oo  “iso -66 


c,  w  .  -  FM 

Cu)  *  [  f  iuu  «  U  (1-  46) 

'looJ-«  'Ju.  rtM 

We  obtain  three  equations  satisfied  by  the  g^(z,u,t)  by  multiplying 

(1.39a)  by  u.  and  integrating  over  du  du  .  These  three  equations, 

1  X  y 

together  with  Eqs.  (1.39b)  written  componentwise,  comprise  the  three 
independent  sets  of  equations  to  be  considered  in  the  next  chapter. 


E  mode 

X 


•3E.  _  _  X 

~  C.  -Jt 


I 

I 

r 


1.  Introduction 

In  this  section  we  present  the  complete  normal  mode  decompo¬ 
sition  for  transverse  disturbances  of  fixed  frequency.  This  speciali¬ 
zation  is  made  in  order  to  illustrate  the  main  features  of  the  method 
which  is  actually  quite  general.  It  will  be  seen  that  a  normal  mode 
decomposition  is  feasible  provided  we  do  not  restrict  ourselves  to 

"functions"  in  the  ordinary  sense  of  the  word,  but  allow  generalized 
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functions  or  "distributions"  as  well.  Except  for  this,  the  normal 
modes  of  the  Vlasov  equation  are  similar  in  character  to  the  normal 
modes  of  any  nonsingular  problem.  They  are  orthogonal  to  the  solu¬ 
tions  of  the  corresponding  adjoint  equation,  and  they  form  a  complete 
set  in  that  any  well-behaved  function  may  be  expanded  in  terms  of  them 
with  the  expansion  coefficients  determined  through  the  orthogonality 
relations. 

2.  Matrix  Notation,  the  Adjoint  Equation  and  Orthogonality 


-19- 


and  take  the  Fourier  time  transform  of  Eqs.  (1.47x).  I.e.  we  seek 
solutions  whose  time  dependence  is  of  the  form  ~  e”^^*  through¬ 
out.  Dropping  subscripts  we  have 


v^3(»,u)= 

f  B(a) 

=.  ^  E^z)  -  J  5fa^tt)ciu 


It  is  convenient  to  combine  the  coupled  Vlasov  and  Maxwell 
equations  into  a  symbolic  matrix  equation  for  the  "state  function" 

^  ,  where  is  written  in  the  form  of  an  array  containing  the  three 


field  quantities. 


g(M)\ 


(2.2) 


(the  z-dependence  is  suppressed) 


Equations  (2.1),  when  combined, 
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become,  in  thi«  notation 


(the  operator  V  is  defined  by 


The  scalar  product  is  defined  in  a  natural  manner.  I.e.  if 


then 


In  this  notation,  the  adjoint  equation  for  the  row  matrix 
ponding  to  ^  becomes 

tfe  = 

7)2  ^  ^ 

with  the  adjoint  operator  defined  by  the  requirement 

By  direct  computation  we  find  that  is  given  by 


H 


t 


o  o 
kii? 

\  C  C 


(2.4) 


(2.5) 


(2.6) 


(2.7) 

corres- 


(2.8) 


(2.9) 


(2.10) 


-21- 

where  the  operator  Jp  is  defined  by 

00 

r  =  f 

‘'F  J-oe 

The  reason  that  we  write  the  adjoint  equation  is  that  solutions  of 
Eq.  (2.4)  and  the  adjoint  equation  (2.8)  are  complementary  in  the 
following  sense:  If  we  seek  solutions  to  Eqs.  (2.4)  and 

(2.8)  with  spatial  dependence  -w  g  ^  then  B/dz  =  i  u/v  and 


(2.11) 


HT„  = 


(2.12) 


(2.13) 


Using  (2.9)  and  the  fact  that  p  is  self-adjoint,  we  obtain 

U.14, 

or 

(T„,,fty)  =  0  V^V  (2.15) 

This  is  the  basic  orthogonality  relation. 

The  orthogonality  relation  will  turn  out  to  be  of  practical  value 
in  computing  expansion  coefficients.  That  is,  the  composite  quantities 
^  are  to  be  thought  of  as  the  normal  modes  of  the  plasma,  and  will 
serve  as  the  basis  for  normal  mode  expansions.  The  orthogonality 


relation  implies  that 


1 
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('f !  f  Ntv)^cv-v') 

V/  J2.16) 

(it  will  be  necessary  to  make  precise  the  meaning  of  the  symbol). 
Then  if 

=  T  A(v)T„ 

V  <2.  17) 

(the  sum  may  run  over  continuous  as  well  as  a  discrete  range  of 
values  of  the  variable!/),  we  inunediately  obtain 

^  ife)  (2.18) 

In  the  following  sections  we  shall  discuss  the  spectrum  of  allowed 
values  of  the  wave  velocity  p  ,  and  obtain  explicit  expressions  for  the 
normal  modes  and  the  normalization  function  N  (V).  Finally,  we 
shall  prove  that  the  form  a  complete  set- -that  an  expansion  of  the 
kind  considered  in  Eq.  (2.  17)  is,  in  general,  possible. 


2.  The  Form  of  the  Normal  Modes 
Let  us  choose: 


ATtiei^c 

OJ 


(2.  19) 


This  choice  is  quite  arbitrary  inasmuch  as  Eq.  (2.  12)  is  linear  and 
homogeneous  in  the  three  field  quantities,  and  is  made  for  convenience.’)' 

The  three  rows  of  the  matrix  equation  (2.  12)  then  become,  after  rearranging: 


It  is  simple  to  see  that  there  can  be  no  non-trivial  solutions  with  By  =  0 
(cf.  eqs.  (2.20)). 
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i^Tziev^ 

60 


(2.20a) 


(2.20b) 


.00 

J  qjXi)(L.  (2.20c) 

-CO 

(W^  =  (4Trne^/m)^^^  is  the  plasma  frequency.)  For  a  given  value  of  u 

£q.  (2.20b)  determines  g^(u)  .  Eq.  (2.20c)  then  serves  as  a  sub- 

oidiary  condition  restricting  the  possible  allowed  values  of  u  .  The 
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solutions  are  similar  in  form  to  those  exhibited  by  Case  for  the 
longitudinal  electric  field  problem.  It  is  convenient  to  classify  them 
in  four  groups: 

Class  1. 


These  are  solutions  for  complex  v  .  We  have 


q^(u;= 

"'’i  '  (j^  ‘  M-Ij 


(2.?^ 


where,  as  required  by  Eq.  (2.20c),  v.^  is  one  of  the  roots  of  the 


characteristic  equation 


{Z.2Z) 


It  is  not  difficult  to  show  that  there  are  only  a  finite  number  of  such 
roots.  Note  that  the  characteristic  function  A(i/)  is  analytic  in  the 
complex  iz-plane  cut  along  these  parts  of  the  real  axis  where  F  (i^)  ^  0  . 
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Class  2  a. 

For  these  solutions  V  is  real  and  F  (v)  ^  0  .  Then  the  solu¬ 
tion  represented  by  (2.21)  is  indefinite  until  we  give  a  prescription  for 
treating  the  singularity  in  the  denonainator.  The  general  solution  of 
(2.20b)  for  real  u  is 

=  +X(V)^(U.-V) 

Co*  W.-V 

The  P  indicates  that  in  integrals  involving  g^(u),  the  Cauchy  prin¬ 
cipal  value  is  implied.  X(i/)  is  arbitrary  and  is  chosen  so  that  (2.20c) 
is  satisfied.  This  condition  becomes: 


A*^  (j^)  and  A  (i/)  are  the  boundary  values  of  A(i/)  as  i/  approaches 
the  cut  from  above  and  below  respectively. 

Clas  s  2  b. 

Here  we  have  i/  real  and  F(i/)  =  0.  The  results  are  essentially 
the  same  as  above:  g^  (u)  is  given  by  (2. 23)  and  X(t')  by  (2.  24).  In 
this  case  (v)  and  A  (j/)  take  on  the  same  value,  since  the  integral 
contained  in  the  definition  of  A(i/)  is  continuous  across  the  real  axis 
when  F(j/)=  0  .  Also,  the  principal  value  sign  is  not  now  necessary, 
but  it  is  carried  along  as  a  reminder  that  when  integrating  with  respect 


to  V  we  are  to  omit  points  of  class  2  c. 
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Class  2  c. 

For  these  solutions  v  is  real,  F  {u)  =  0  and  \{u)  =  0  as 
well.  Just  as  in  class  1,  there  are  only  a  finite  number  of  u, 
satisfying  these  conditions.  The  solution  for  g  (u)  is 


(2-25) 

If  F  (u)  vanishes  at  least  linearly  as  u  g  (u)  is  a  per- 

H 

fectly  well-behaved  function.  We  shall  assume  this  to  be  the  case. 

At  any  rate,  the  (possibly  improper)  integral  /  g^  du  is  well 
defined. 

To  recapitulate:  the  fundamental  equation  (2.12)  possesses 
solutions 


t  - 


\ 


\6,/ 


(2.26) 


of  two  basically  different  types.  There  is  a  discrete  set  of  values 

/^j)  =  o 


I/.,  V.  with  either 
J  1 


a)  complex, 


or  b)  i/.  real,  /^(v.)  -  P  )  =  0 

The  respective  g^  (u)  and  g^  (u)  are  well-behaved  functions.  The 
J  i 

remaining  (continuum)  real  values  of  v  are  associated  with  solutions 
for  which  the  g|^(u)  are  singular  and  are  interpreted  as  distributions 
in  the  Schwarzian  sense. 


-26- 


3.  Solutions  of  the  Adjoint  Equation 

For  every  allowed  eigenvalue  of  the  fundamental  equation  (2.12), 
there  ia  also  a  solution  of  the  adjoint  equation  (2.13).  The  decom¬ 
position  into  classes  is  essentially  the  same. 

Class  1  :  v  complex 

We  take,  for  convenience 


t  -  nev.c 
miu) 

Then  the  three  columns  of  the  matrix  equation  (2.13)  read 


(2.27) 


mitd 


CU-Vy)  9^.(U)  =S 


U)*- 


v- 


The  Solution  of  (2.28b)  is 


J*',-  W*  j  M' 


1- 


which  when  substituted  into  (2.  28  c)  yields  the  condition 


A('')  =  0 

as  before. 

Class  2  a.  v  real  and  F  {v)  ^  0 


(2.  28a) 


(2.  28b) 


(2.28  c) 


(2.29) 


(2.30) 


Again  we  choose 
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-t  neuc*‘ 

(2-27> 


which  implies 


mi  u) 


{2.28a) 


Equations  (2.28b,c)  are  also  still  valid,  but  the  solution  for  this  cas'* 


is 


*  XMS(u-v)  (2.31) 


with 

Xiv) 

F(v) 

Class  2b.  1/ real,  T  (v)  =0,  but  X(i/)  ^  0 


Here  we  must  take 


(2.32) 


A  nonzero  choice  for  E  must  be  avoided. 


(2.  33) 

It  would  lead  to  solutions 


V 

of  the  same  form  as  those  of  class  2  a.  But  this  is  impossible,  since 
Eq.  (2.  32)  cannot  be  satisfied.  The  matrix  equation  (2.13)  becomes: 


(2.  34a) 

(2.  34b) 


FM^I('u)cLu 


=  0 


(2.34  c) 
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The  solution  of  Eq.  (2.34b)  is 

which  when  substituted  in  Eq.  (2.  34  c)  yields  an  identity. 
Class  2c.  real,  F(i/^)  s  0  and  X(j/^)  s  0 
As  before,  we  take 


which  implies 


•nico 


(2.35) 


(2.27) 


gt  _  y\6V^c 

ho  i  u> 


(2.28a) 


J  F (■w)  <Au  s.  c  *■-  V.  *■ 

-so  * 


(2.28b) 


(2.28  c) 


The  solution  for  g  (u)  is 
^i 


0). 


(2.  36) 


4.  The  Normalization  Coefficients  N(i/) 

The  general  orthogonality  relation  (2.15)  may  be  verified  by 
direct  calculation,  using  the  solutions  obtained  in  the  preceding  sec¬ 
tions.  It  remains  to  define  precisely  and  determine  the  normalization 
coefficients  N(|/)  referred  to  in  Eq.  (2.16).  It  is  simplest  to  consider 


the  discrete  spectrum  first. 


-29- 


Suppose  we  have  a  .  Then  the  orthogonality  relation  implies 


and  since 


(2.15') 


^  r\eV;C  MTciei^.’’  ^  nei^‘c  /#tii€t/c 

r»i  i  </j  u>  Witt;  to 

tt)‘  J  L  tt)^  J  J-eo  J 


which  becomes  after  a  little  algebra, 


i  AM 

40*  J 


Vssl'- 


Thus  we  say 


=  NlvpS,_v. 


with 


'V*1J 


(2.37) 


(2.38) 


(2.39) 


I/, 

«l 


(2.40) 

is  the  Krone  eke  r  6 -symbol.  Since  i'.  is,  by  assumption. 


a  simple  zero  of  the  characteristic  equation,  N(i/j)  will  not  vanish 


(cf.  (2.18)). 
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The  result  is  essentially  the  same  when  we  consider  a  . 


(2.39') 


N(v)  = 


60*-  ‘ 


(2.41) 


T***  1 
60*  ‘  ^ 


iAM 


V-W^-lfe  p 


When  i>  belongs  to  the  continuum,  the  situation  is  somewhat 
different.  The  corresponding  normalization  integrals  are  undefined. 
We  write,  however, 


(2.42) 


where  the  (Dirac)  6 -symbol  is  nothing  more  or  less  than  an  abbrevia¬ 


tion  for  the  following  statement;  if 

-b 

*  I  -oo4a<;bioo  (2.43) 


then 


/  r  nTr\  rA(v)N(v)  a<v<b 

^  0  Olberwi^e 

It  is  a  simple  (and  familiar)  matter  to  calculate  N(  {/)  vdien  v  belongs 
to  class  2b.  Then  ^  is  of  the  form 
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6v(’*)A(v')clv' 

a 

f**  E^,  A(v)<<v' 

•'<1  ^ 

E>^.  A(v';<iy' 


(2.45) 


But  since  E  =  B  =  0  we  have 
1/  1/ 


_  r  5^t^u)uQ(u)<iu  =  f  gJ^u^utiuT  %,(u)A(v')cIu' 

-CO  -M  -^A 

(2. 

=  r^^(u-v)4iA|  [X(u)A(u)  a<«<t') 

“00  L  -la  0  otherwise  / 


The  first  term  in  the  braces  does  not  contribute  to  the  final  result 
since  F{i/)  =  0  .  Thus 


(%Nt) 


vA(i/)A(v)  «<v<b 

0  otherwise  (2.47) 


which  implies 


N(V)  =•  I^A(V)  V  m  class  Zk.  (2.48) 

The  preceding  result  is  little  more  than  a  partly  rigorous  justification 
of  the  often  used  symbolic  relation 


I  f(-u.)  SlU~V)  SiM-vOd-U  -  flv)S(V-l/‘) 

-  Pi 

but  it  is  important  to  make  precise  the  meaning  of  the  normalization 
for  the  continuum  functions.  This  requires  that  we  consider  not  the 
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functions  themselves  but  "wave  packets"- -superpositions  of  the  ^ 
where  u  varies  over  a  continuous  range.  Equations  (2.43)  and 
(2.44)  follow  immediately.  It  would  be  tempting  to  conclude  (incor¬ 
rectly)  from  (2.42)  that 

Hlv)  ^  a<v<y  (2.50) 

but  in  fact  this  is  just  Eq.  (2.44)  with  the  orders  of  integration  inter¬ 
changed.  This  point  is  emphasized  because  when  i/  belongs  to  class 
2  a,  we  come  upon  one  of  the  rare  instances  that  care  must  be  taken 
to  perform  the  integrations  in  the  correct  order. 

In  calculating  N(i/)  for  j/  in  class  2a,  we  note  that  there  is 
no  objection  to  letting  a  =-  b  =  in  (2.43),  just  so  long  as  we 
make  sure  to  choose  A(i/')  so  that  all  the  integrals  exist.  Assuming 
this  to  be  done,  we  have, 

A(V)N(W  = 

+  A(W5(w-l/)jwelu  I  r^^V'^P-;^^^+X(v')^(u-v')jA(v')flJv' 
+  r”4Tciev'M(v‘)clv*  n  ev*’c  4nte  v'c  A(v‘)dv' 

J  60  mito  60 

-0) 

=  (ypfr»  pf“  P 

\WV  'Lao  “M-  - 


r- 

^-00 


V*’ A(v')olv' 

V'  -M. 


-h 
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(A)*'  IX-V  w*  J_^  V‘  -1/ 


+v 


F(v) 


£i)  *  r  ^ 

A(v)  +  (vv'^+v'v*-)  A(i/')dv' 

to*  I 

-00 


(2.51) 


The  first  term, containing  the  iterated  principal  value  integrations, 


is  the  only  one  that  presents  any  difficulty.  Here  we  use  the  Poincare- 
Bertrand  formula^^  to  obtain 

^"uF(u)c(u 


J-oo 


Mv‘)<^v' 


(2.52) 


“  (^  ^*)  ^  {^*)  v'^MvOdv'  pj  , 


U  Ru)  du 


Jtt-v)(u-v') 


The  integration  over  u  may  now  be  performed  using  partial  fractions 
and  the  definition  of  X(l/)  .  The  result  exactly  cancels  all  the  remain¬ 
ing  integrals  in  Eq.  (2.51)  leaving 


mv)  = 


X^v) 


\v) 


Fiv) 


(2.53) 


V  in  class  2  a. 

5.  The  Normal  Mode  Expansion  and  the  Full  Range  Completeness  Theorem 
Here  we  prove  that  the  normal  modes  form  a  complete  set, 
in  the  sense  that  an  arbitrary  state  function  ^  may  be  expanded  in  terms 
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r 

I 

I 


oi  them. 


1  j  *  J-eo 

Once  we  know  that  the  expansion  is  possible,  then  the  expansion 
coefficients  are  determined  using  the  orthogonality  relation  (2.15) 
and  the  known  normalization  functions,  viz. 


(2.54) 


with  N(i/.)  given  by  (2. 41) 


’  N(v.)^  J*'  N(v)  j_  60’  ^  'iflo  vniui  L  ^  -!^2. 

with  NQ/.)  given  by  (2, 40).  When  i/  is  in  class  2  a, 

J 


l<v)  9(V)  t  Flv)  tt'pl 

1t»[^V)V’(v;  +A‘(v) 


(2.57) 


and  this  expression  is  also  correct  for  i/  in  class  2  b,  since  in  this 


Ai’') ' 


(2.58) 


which  is  just  Eq.  (2.57)  with  F  (j/)  =  0  . 
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The  full  range  completeness  theorem  can  be  simply  stated  as 

follows; 

THEOREM  I. 

Given  a  state  function  ’S'  with  g(u)  everywhere  well-behaved, 

/  g(u)du  <  but  otherwise  arbitrary,  then  an  expansion  of  the 
form  presented  by  Eq.  (2.  54)  exists,  and  the  expansion  coefficients 
a^,  a^  and  A(t;)  are  unique. 

In  the  proof  "well-behaved"  will  be  taken  to  mean  "satisfying 
a  Hdlder  condition".  That  is  there  exists  a  pair  of  numbers  M  and 
y  (y  >  0)  such  that  for  any  u^ 

I  g(u)-  0(llo)|  <  (2.59) 


Actually,  the  theorem  holds  for  a  much  wider  class  of  functions  g{u). 
For  example,  g  (u)  may  itself  be  a  Schwarzian  distribution.  In  the 
same  spirit  we  shall  assume  that  F  (u)  also  satisfies  a  Holder  condition. 

The  proof  of  the  theorem  involves  an  actual  construction  of  the 
expansion  coefficients.  That  is,,  we  show  that  (2.  54)  possesses  a  solu¬ 
tion  by  solving  it.  Write 


i  •  j  > 


(2.60) 


with  given  by  (2.55,56) 


We  note  that  the  expansion  of  ’i' '  in  terms  of  the  continuum  functions 
alone 

r  ^ 

■  I  A  ,  \  nir  J.. . 

(2.61) 


pOO 

*Y'  «  I  A(v)Tvdv 
'^-00 


is  equivalent  to  the  expansion  (2.54)  of  ♦  in  terms  of  all  the 
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^ '  must  satisfy 


•V  *  ) 


(2.62) 


as  well  as  the  conditions  imposed  on  ^  in  the  statement  of  the  theorem, 
but  is  otherwise  arbitrary.  The  three  components  of  Eq.  (2.61)  are 


g',u)=  +A(»)ACU) 


(2.63a) 


E'  = 


tjTcte 

Ui 


f  V^A(-^)dv 

^  J.eo 


(2.63b) 


B  = 


f*®., 

“L 


=  I  vAMdv 


(2.  63c) 


Equation  (2.63a)  is  a  singular  integral  equation  of  a  type  treated 

extensively  by  several  Russian  mathematicians.  The  solution  closely 
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follows  Muskhelishvili.  We  shall  make  use  of  the  following  properties 
of  singular  integrals  of  the  Cauchy  type: 

Suppose  M  (x)  satisfies  a  Holder  condition  and  dx'  <  «> 

then 


a) 


M(x*jcU* 

X'-i 


is  an  analytic  function  in  the  complex  z -plane  cut  along  the 
real  axis  form  a  to  b,  vanishing  as  z  at  least  as 

fast  as  l/z  . 

b)  7>\^(x)  and  ')71*(x),  the  boundary  values  of  171  (z)  as  z 
approaches  the  cut  from  above  and  below  respectively, 
each  satisfy  a  Holder  condition. 
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c) 


7n^(x)  -'rrcit)»  M(x) 
Tci  [yrfu)  -t-yiTM) «  ?J 


{2.64a) 

(2.64  b) 


(the  Plemelj  formulae) 


A 


vanishes  as  z  -►  ®® 


The  converse  also  holds:  if  yn  (z)  satisfies  a)  andb),  then  it  may  be 
represented  as  a  Cauchy  integral  of  the  form 


x'-i 


(2.65) 


M(x)  will,  of  course,  be  given  by  the  first  Plemelj  formula. 

We  combine  {2.63a-c)  into  a  single  equation  for  A{u)> 


ug'(u)  =  - 


+X(u)x.A(u) 


Now  define 


xt  0'(u)dlu 

u  -a~ 


(2.66) 


2 


V  A(v)dv 


V  -2 


(2.67) 


(2.68) 


(2.69) 
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(z)  and  ^ (z)  are  analytic  in  the  cut  plane  and  vanish  at  infinity. 
Their  boundary  values  (u)  and  Q  —  (li)  satisfy  a  Holder  condi¬ 
tion.  If  a  solution  A(j/)  to  (Z.6b)  exists,  then  ^(z)  will  have  those 
same  properties.  We  shall  show  that  such  an  A(i/)  exists  by  con¬ 
structing  the  appropriate  di^)  a^nd  showing  that  it  does  in  fact  possess 
the  properties  we  require  it  to  have.  We  would  then  find  A(i/)  by 
using  the  Plemelj  formula 

AM  =  i 

if  we  didn't  already  have  the  orthogonality  relations.  In  terms  of  the 
boundary  values  of  the  functions  ,  'J'  and  d  •  Eq.  (2.  66)  becomes, 
after  some  rearranging  and  dropping  primes 

^  (2.71) 

-  aV)] y'M -(c"-u*)a'(w 

Now  consider  the  function 

Jh)=  (2.72) 

J  (z)  is  analytic  in  the  cut  plane,  since  f  and  5  are,  and  d  is 
assumed  to  be.  But  according  to  (2.  71),  J  (z)  is  continuous  across 
the  cut.  Therefore,  J (z)  is  an  entire  function,  and  may  be  found 
simply  by  observing  its  behavior  for  large  values  of  |  z|  We  note 


that 
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O 


1*1 -*^00 


(2.73) 


F('^<lu 


“/*  |«/ 


■»  - 


airi 


(2.74) 


-  ZKi  + 


+  -L-  f _  — L  JiL  fF+Af(\^l, rVA(v)<iu 

Zni)  V-2  a,ili  Wiel^^‘=''V  iiU  -— -- 

-M  *• 


Itl-ttO  3Xi 'flu# 


(2.75) 


(f*|8) 


(the  fact  that  /  y  A{i/)  di/  <  «  will  be  verified  later).  Thus 


J<») 


«  ■■  ii^r  ^  ^ 


|2l— ^eo  SrKl  <0[ie 

We  conclude,  from  Liouville's  theorem 


(2.76) 


J'*>= 

which  gives  the  solution 


(2.77) 


at*) 


c‘-2^  -i1Ci^‘z’T(*) 


(2.78) 


Now  we  must  examine  the  function  (\(z)  and  see  if  it  has  the  proper¬ 
ties  demanded  of  it.  We  note  that  C^(z)  vanishes  as  |  z  |  -•  ®*  and 
is  analytic  in  the  cut  plane  except  possible  at  the  zeroes  of  the  denom¬ 
inator,  where  ^(z)  could  have  simple  poles.  But  these  points  are 


just  the  V.,  and  when  p  =  v.  ,  the  numerator  may  be  shown  to 
J  J 


vanish,  viz. 
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[  numerator]  X 1/*  2Ki  ^ 
J  (*) 


(2.79) 


Thus  the  singularity  is  removable,  and  OiM  analytic  in  the  cut 
plane.  We  need  look  only  at  its  boundary  values.  As  z  approaches 
the  cut  from  {  above  ■»  denominator  becomes 

[denominator]  =  c"-- V*  -  V^Xl[7*(v)  *  ^  (2-  80) 

The  boundary  values  are  well-defined  except  possibly  where  the  denom¬ 
inator  vanishes,  so  the  points  in  question  are  just  the  t/^.  And,  at 
these  points,  the  numerator  also  vanishes,  since  by  direct  calculation 


[  nvimerator]  X  ll'aiTCi  “  0  (2.81) 

Finally,  we  must  show  that  the  various  integrals  involving  A{i/)  actually 
exist.  By  direct  calculation,  we  find 


CO  R 


I 


9IC{CC  * 

This  implies  (cf.  condition  a),  page  36)  that  /  A (i/)  dp  <  . 

(It  is,  in  fact,  given  by  (i;B/4ir  ie  c.)  Moreover,  by  condition  d),  the 
Cauchy  integral 


(2.82) 


a(.)  =/" 

•00 


(2.83) 
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converges,  and  another  straightforward  calculation  yields 


.00 

au)  =  ^au)  VA(v)o(v 


_ ^  I  j_ 

-to  !*!-♦«)  inciec  XKi  z- 

which  implies  that  /  A(^')  di/  converges  and  its  value  is 

wE/ 4 Trie  .  Finally,  we  note  that  this  means  that  the  Cauchy 
integral  /  y  A  (u)/ (i/ -  z)  dj/  converges.  Hence  0{{z),  the 
solution  of  (2.71)  does  indeed  possess  the  properties  we  assumed 


it  had.  This  completes  the  proof. 

For  future  reference  we  note  that  by  virtue  of  the  Plemelj 
formulas  N{ir)  may  be  expressed  in  terms  of  the  boundary  values 
of  the  analytic  function  A  (z) 

^  _ nv) _ 

N  M  v[x(v) +irt  V*  fiv^[\{V)  -xi  ^!p(vJI 


_  F(V)  _ 

V  A'  (v)  3X1  *  /\*(v)A'‘M 

=  7'(^) 

where 


ixi  s'A  (*) 

6.  Specialization  to  the  Boltzmann  Distribution 

We  consider  now  the  particular  case  where  (u)  is  the 
Boltzmann  distribution  for  electrons  at  temperature  T  =  @/k. 


(2.84) 


(2.85) 


(2.86) 
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-1—1 

(zTceJ 


(2.87) 


Then 


/  iM  i/l 


(2.88) 


Some  of  the  results  of  this  section  will  depend  only  on  the  general 
shape  of  F  (u)  and  not  on  its  exact  functional  form. 

It  is  possible  to  ascertain  the  number  of  elements  in  the  dis¬ 
crete  spectrum;  i.  e.  the  number  of  roots  of 


AW  =  c  -  w>*  J  -i(-*  -  0 

-oO 


(2.89) 


First  note  that  since  A(z)  is  a  real  function  of  z,  the  roots  occur 
in  complex  conjugate  pairs.  But  since  F(u)  is  an  even  function,  A(z) 
is  also  even  and  -  z^  is  a  root  along  with  z^.  Thus  in  general,  the 
roots  occur  in  groups  of  four;  z^  and  the  upper  half  plane, 

symmetric  about  the  imaginary  axis,  and  -z^  and  z^*^  in  the  lower 
half  plane.  We  shall  show  that  there  is  at  most  one  root  in  the  upper 
half  plane.  This  root,  when  it  occurs,  must  then  lie  on  the  imaginary 
axis. 


We  seek  to  find  the  number  of  zeros  of  A  (z)  in  the  upper  half 
plane,  using  the  argument  principle.  Consider  the  contour  C  in  the 
upper  half  z-plane.  (See  Fig.  1.)  Let  us  follow  the  behavior  of  A(z) 
as  z  traverses  the  contour  C.  The  change  of  argument  of  A(z)  will  be 
equal  to  2ir  X  (number  of  zeroes  of  A(z)  inside  C).  The  semicircle 
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1-2-3  is  large  enough  so  that  A(z)  assumes  its  asymptotic  behavior 
along  it. 


For  z  on  the  large  semi-circle 

A(*)a  -  (i- (2- 90) 

and  when  z  -  x  4  ic  is  just  above  the  real  axis,  we  have 

A(*)=  A"(«)= 


Mo)  =  C'’ 


(2.92) 


Since  F(x)  is  positive  and  nonvanishing,  when  z  is  on  the  upper  lip 
of  the  real  axis.  ImA(z)  is  for  {"'f 

vanishes  only  when  z  -  0  .  Thus,  as  z  traverses  the  path  3-4-1, 
A(z)  crosses  the  real  axis  only  once.  There  are,  then,  only  two 


possibilities. 
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} 


1)  U)  >  U)  The  image  of  the  contour  C  is  shown  in  Fig.  2. 

P 

A  arg  A(z)  =  0  and  there  are  no  zeroes  of  A(z)  inside  the 
contour  C;  hence  no  zeroes  in  any  finite  part  of  the  upper  half  plane. 


A  argA(z)=  2  77  and  there  is  one  zero  in  the  upper  half  plane. 


Thus,  for  Boltzmann  F(u),  the  set  of  allowed  values  of  if 
includes  the  entire  real  axis  (class  2  a)  and  sometimes  the  two  dis¬ 
crete  roots  +  it'g  •  (class  2  b) 


t 
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We  may  also  obtain  an  estimate  of  the  value  of  the  discrete 


root  ii/^.  We  have,  from  Eq.  (2. 89)  using  the  explicit  form  of  F  (u) 


(2.93) 

Assume,  for  the  moment,  that  | 

z  1  is  large  enough  so  that  the 

asymptotic  expression 

-  1 

(2.94) 

X-/^2 
-00  Hzo 

is  valid.  Then 

(2.95) 

or 

C*' 

CO 

_ _ _ 

(2.96) 

z  = 

1  -  HVco*-  , 

This  result  is  correct  so  long  as  (2.94)  holds,  which  amounts  to 

I© 


©  >1  *  5 


(2.97) 


or 


^  c| 


* 

which,  for  T  »  6000°K,  becomes 


(2.98) 


CO  >  4  X  10' 


(2.99) 
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Below  this  frequency,  we  may  try  to  approximate  the  result  by  assum¬ 
ing  that  I  (m/2  6)^^^  z  |  turns  out  to  be  small.  Then  we  may  use 
the  power  series  expansion 


U® 


(2.100) 


u® 

To  lowest  order,  this  yields 


V  = 
o 


(0 


(2.101) 

and  an  iteration 


lives 


(2.102) 


We  note  that  except  at  extremely  high  temperatures  or  very  low  fre¬ 
quencies,  the  location  of  the  two  modes  with  imaginary  wave  number 
is  correctly  given  by  the  zero  temperature  treatment. 

Finally,  we  may  obtain  an  approximate  expression  for  the 

normalization  coefficient  N  (  +  i  )  defined  by 

—  o 


N  (±iv^  =  -  %  «  am 


*  a 


CO 


dv 


By  direct  calculation,  for  Boltzmann  F(u),  this  becomes 


(2.103) 
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N(±iv„)=  -  (2- 


104) 


Except  at  very  low  frequencies,  to  first  order  in  2  0/mc  we  get 


Nr±iH)=  tzl 


60^ 


-  Tl  —  ^  -^1  ^  me  V 

*•  [_'  4a)*  /  a;L»  ,  \Vi 


(2.105) 


Chapter  III.  BOUNDARY  VALUE  PROBLEMS 
FOR  TRANSVERSE  WAVES 


1  ■  Introduction 

In  this  section  we  shall  utilize  the  normal  modes  of  the  Vlasov 
equation  in  the  solution  of  some  typical  boundary  value  problems.  The 
method  is  straightforward:  we  write  the  state  function  of  the  medium 
under  consideration  as  a  superposition  of  normal  modes,  and  apply  the 
boundary  conditions  to  determine  the  e:q>ansion  coefficients.  In  certain 
cases  (the  degenerate  problem)  the  expansion  to  be  perfornr^d  will  utilize 
the  full  set  of  normal  modes,  and  the  orthogonality  relations  may  be  used 
to  determine  the  coefficients.  In  any  event,  we  will  arrive  at  an  integral 
equation  for  the  expansion  coefficients,  part  of  which  (the  so-called 
"dominant"  part),  by  virtue  of  the  singular  nature  of  the  g^  (u),  will  con¬ 
tain  singular  integrals. 

We  shall  consider  the  reflection-transmission  problem.  An  electro¬ 
magnetic  wave  of  frequency  u  is  incident  on  a  plasma.  What  is  the  strength 
of  the  reflected  wave  ?  It  will  be  convenient  to  ignore  the  fact  that  the 
actual  plasma  oscillations  are  "forced"  by  the  impinging  wave.  This  avoids 
the  problem  presented  by  the  presence  of  an  external  force  field  whose 
effect  has  not  been  taken  into  account  in  the  linearization  of  the  basic  equa¬ 
tions.  The  electromagnetic  fields  throughout  the  plasma,  in  particular  at 
the  boundary,  are  to  be  viewed  as  "small"  in  the  sense  of  the  original  linear¬ 
ization  (Chapter  I).  We  simply  construct  a  solution  representing  a  plasma 
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configuration  consistent  with  the  presence  of  an  incident  and  reflected 
wave,  and  only  afterwards  do  we  interpret  the  plasma  motions  as  the 
effect  and  the  impinging  wave  as  the  cause. 

2.  Boundary  Conditions 

The  boundary  conditions  for  the  electromagnetic  fields  are  simple. 
Since  all  disturbances  are  small,  there  can  be  no  surface  changes  or 
currents  and  the  fields  must  be  continuous. 

£  (2  =  0*)=  E(b  =  0-) 

(3.1) 

B  f2=  0+)  -  £(2  -0-) 

We  must  specify,  in  addition,  the  behavior  of  the  electron  distri¬ 
bution  function  at  the  boundary.  It  is  simplest  to  assume,  with  Landau*,* 
that  electrons  experience  specular  reflection  back  into  the  plasma  at  the 
boundary.  This  automatically  implies  that  at  the  boundary  the  distribution 
function  is  symmetric  in  the  normal  component  of  the  velocity 

or 


q(OjU)  =g(o,-u; 


(3.3) 
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This  is  the  degenerate  case,  and  the  solution  may  be  written  in  terms 
of  the  normalization  integrals  using  the  orthogonality  relations. 

The  opposite  extreme,  one  might  term  the  "diffuse  absorber.  " 
In  this  model,  the  electrons  are  returned  to  the  plasma  at  the  boundary 
in  a  manner  essentially  independent  of  the  velocity  distribution  of  the 
electrons  which  leave.  This  implies; 


fl(lA,2=o)  =  (a  given  function) 

*  •14>0 

This  boundary  condition  leads  to  a  singular  integral  equation  which  is 
equivalent  to  the  inhomogeneous  Hilbert  problem. 

The  intermediate  case,  too  complicated  to  be  considered  here, 
allows  the  distribution  function  of  the  incoming  electrons  to  depend  on  the 
outgoing  electron  distribution  in  an  arbitrary  (but  linear)  manner.  This 
leads  to  a  singular  integral  equation  for  the  e3q>ansion  coefficients  which 
contains  a  Fredholm  term  as  well. 

3.  The  Electromagnetic  Fields  at  the  Boundary 

We  consider  a  plane  polarized  electromagnetic  wave  impinging 
from  the  left  on  the  surface  of  a  plasma  (interface  at  z  =  0).  In  Chapter  I 
it  was  pointed  out  that  the  transmission  and  reflection  coefficients  for  an 
ordinary  dielectric  medium  may  be  written  in  terms  of  the  index  of  refrac¬ 
tion  of  the  medium,  (cf.  Eq.  (1. 6). )  Of  course  in  this  general  treatment, 
the  plasma  is  not  simply  a  dielectric  medium.  For  a  given  frequency,  waves 
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of  all  wave  numbers  can  propagate.  Nevertheless,  it  is  an  easy 
matter  to  define  a  quantity  which  plays  the  same  role  as  the  index  of 
refraction. 

The  incident  wave  is  of  the  form 


(3.5) 


The  reflected  wave  is  of  the  form 


(3.6) 


Hence,  by  (3. 1)  the  electric  and  magnetic  fields  at  the  boundary  are 
given  by 


£■^(2-0) 


B,  (2-o)=6„e-‘<(£;-£.)e-'“^ 


(3.7) 


Thus 


£r  ^  LzA 

£.  I -i- n  (3.8) 

where  n  is  the  ratio  of  the  magnetic  and  electric  fields  at  the  boundary 


(3.9) 
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Note  that  this  is  a  reasonable  analogy  to  the  index  of  refraction  of  an 
ordinary  dielectric  medium,  since  if  we  assume  a  disturbance  in  such 
a  medium  to  be  of  the  form 


(3.10) 


then  it  follows  from  Maxwell's  equations  that  n  represents  the  ratio  of 
magnetic  and  electric  fields  throughout  the  medium. 


4.  Reflection  From  a  Plasma  Half  Space 

We  consider  the  state  function  of  the  plasma  to  be  expanded  in 
terms  of  the  normal  modes 

,3. 

V 

The  summation  incorporates  an  integral  over  the  continuum  modes  and 
a  discrete  sum.  The  boundary  condition 


a(o  u)  gfo, -w) 

J  (3.12) 


implies  that 


(3.  13) 
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We  will  deal  here  only  with  the  case  treated  in  Chap.  II,  sec.  6, 

where  F  (u)  is  assumed  to  have  such  a  form  that  there  are  either 
o 

no  discrete  eigenvalues  (when  u  >u^)  or  two  imaginary  eigenvalues 
=  ±  i  (when  u  <«p)-  case 

(3.14) 

for  both  the  discrete  and  continuum  modes,  and  hence  from  (3. 13) 


^  [A(v)-A(-v)]^^(u)^0 

V 

Moreover, 


(3.15) 


5£a(W-A(->')J£'^  =  o 

V 


simply  because  E  is  even  in  v  .  Also,  since  B  is  odd  in  v 

V  V 


(3.16) 


V 

Therefore,  if  we  define 


(3.17) 


C(u)  =Mv)-M-v) 


we  obtain  the  concise  relation 


(3.18) 


(3.  19) 
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Sctv)%  ~ 


which  implies  that  C  (  v)  is  the  coefficient  in  the  expansion  of 


(3.20) 


in  terms  of  the  i|j  Hence,  by  the  orthogonality  relation  (2.  18) 


Civ)  * 


a:) 

N(v) 


(3.21) 

Having  obtained  C  {  v),  in  order  to  obtain  the  desired  coefficient 
A  (  v)  we  make  use  of  the  fact  that  the  system  must  be  well  behaved  as 


(3.22a) 


y)  f\W  -  0  V<0 


(3.22b) 


Condition  a)  insures  that  we  exclude  the  exponentially  growing  modes, 
and  b)  may  be  thought  of  as  a  radiation  condition  (  no  incoming  waves)  or 
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a  regularity  condition  which  results  if  we  assign  to  u  a  small  imaginary 
part.  This  last  point  will  be  discussed  in  section  7  of  this  chapter. 

Combining  (3. 18),  (3.21)  and  (3.22),  we  obtain  the  estpansion 

N(y) 

v<o 

(3.23) 

A(^^o)  *  ^ 

from  which  g  (z,  u),  and  £  (z),  and  B  (z)  may  be  computed.  In  particular, 
since 


E  =  5 


(3.24) 


we  obtain  for  the  index  of  refraction  n 


JL  — 

^  * 


6. 


=  Z 


f  bJe^cIu  ^ 

Jo  Ni-^)  NiMU 


(3.25) 


[ 


[ 


where  the  discrete  term  is  to  be  included  only  when  appropriate  (i.e.  when 

u  <  u  ). 

P 
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Using  (2.85-6),  the  integral  term  in  (3.25)  may  be  written  as 
a  contour  integral. 


Moreover,  in  view  of  (2.40),  the  discrete  term  becomes  simply  the 
residue  of  [  -  2c/A(z)]  at  z  ^  -  i  v^.  This  leads  to  the  compact  expres¬ 
sion 


(3.27) 


In  the  zero  temperature  limit 

F(u)  =  S(u) 


(3.28) 


This  implies 

A(2)  =  C  ”  (  ^  To’O  (3.29) 

and  the  contour  integration  may  be  performed  easily,  to  yield  the 


familiar  result 
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n 


U) 


(3.30) 


The  zero  temperature  result,  which  predicts  a  purely  imaginary 
index  of  refraction  (and  corresponding  complete  reflection)  for  frequencies 
below  the  plasma  frequency,  has  well-known  experimental  consequences. 

It  provides  the  es^lanation  for  the  "plasma  blackout"  which  inhibits 
communication  with  a  reentering  space  vehicle.  It  is  evidenced  also  in 
the  relative  ease  with  which  the  low  frequency  AM  signals  may  be  trans¬ 
mitted  over  large  distances  by  reflection  from  the  ionosphere,  whereas 
transmission  of  the  higher  frequency  FM  signals  is  limited  by  the  curvature 
of  the  earth.  Nevertheless,  at  finite  temperatures,  since  the  integrand  in 
(3.26)  is  explicitly  positive  definite,  there  will  be  a  small  contribution  to 
the  real  part  of  the  index  of  refraction  at  frequencies  below  the  plasma 
frequency,  the  real  expression  (3.26)  serves  as  a  slight  modification  to 
the  zero-temperature  result  n  =  (1  - 

For  purposes  of  computation,  it  is  convenient  to  write  (3.27)  in 
the  form 


I  f  di 

vl  2tiJ 


(3.31) 
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where  U  is  the  velocity  of  light  in  units  of  (2  and 

A''(Z)=  = 


-00 


(3.32) 


Z  (z),  the  "plasma  dispersion  function,  "  is  tabulated. 
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5.  Reflection  From  a  Slab 

We  consider  *’he  reflection  of  an  electromagnetic  wave  of 
frequency  u  incident  on  a  slab  of  thickness  a.  As  in  the  case  of 
a  half  space,  we  write  the  state  function  as  a  superposition  of  normal 
modes  and  apply  the  boundary  conditions.  We  again  obtain: 

C(v)^A(^)-A('^)  =  2Bo 

(3, 

We  must  also  apply  the  condition  of  specular  reflection  to  the  far  side 
of  the  slab. 


g  (a,u)  = 


(3.34) 


This  gives 


D(v)- 


i-a  -iv® 

A(v)e''-AMe  =26«^|^ 


(3.35) 


where  B  is  the  magnetic  field  at  the  surface  »  =  a, 
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Combining  (3.33)  and  (3.35)  we  obtain  the  expansion  coefficient 

...  /,  f  R  ^  R 


(3.36) 


which  we  use  to  obtain  the  electric  fields  at  the  boundaries  in  terms 
of  the  magnetic  fields 


=  oiB^  +jSB, 


(3.37) 


-  $  A(v)e‘“‘*£  -  s  i S ~  B.i S 


®  ^  Nlu)slfi  ^ 


ft  (  \  ^ 

'vN(i))tary^ 


=  -1^4  - 


(3.36) 


The  above  expressions  have  been  simplified  by  the  use  of  the  fact  that 
B  and  N  (  v  )  are  even  and  odd  functions  of  v  respectively.  The 
last  boundary  condition  to  be  applied  follows  from  the  requirement  that 
the  electromagnetic  field  in  the  space  behind  the  plasma  slab  represents 
an  outgoing  monchromatic  wave 


Combining  the  la|t  three  equations,  we  obtain  the  result 


(3,39) 


(3.40) 
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The  singularities  in  the  integrals  defining  o  and  ^  are  to  be  resQlved 
by  ascribing  to  u  a  small  positive  imaginary  part,  as  will  be  discussed 
later. 

Therefore,  with  respect  to  a  normally  incident  plane  wave,  the 
slab  behaves  like  a  half  space  with  index  of  refraction  given  by 


The  coefficients  a.  p  may  be  simplified  via  contour  integration 

_ Lu)a. 

-to 

(3.42) 


«  I 

Z 


Do 


(-1)" 


=  yP-Z. 


(3.43) 


For  the  zero  temperature  model  [  A(z)  s  c*  -  «  z*  ]  ,  we  obtain 


t 


(3.44) 


which  is  in  agreement  with  the  result  obtained  by  elementary  methods. 


-61- 


6.  Reflection  From  a  Half-Space:  The  Diffuse  Absorber 

The  problem  of  the  diffuse  absorber  is  of  ixiterest  because  it 
illustrates  a  situation  where  the  eigenfunction  expansion  must  be  per¬ 
formed  without  the  aid  of  an  orthoganality  relation. 

We  assutpe  that 


•»  a 

•21I& 

ze 


(3.45) 


which  implies 


0(0^  it)  =  O 


^>0 


(3.46) 


Then  with  u  >  u  so  that  there  are  no  discrete  v.  •  the  eigenfunction 

P  J 

expansion  becomes 


OO 


1A,>0 


(3.47) 

Thus,  we  are  presented  with  a  "half-range"  expansion--the  expansion 
of  (o  <  u  <  co)  in  terms  of  the  ij,  (o  <  y  <  oo).  The  resulting 
integral  equation  is 


0  =  -  R«)  ^  A(W  « >  0 


(3.48) 


subject  to 
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J  V  Alv)civ  * 

■^0 

A3  6^ 
UtiZC 

f  v’-A(v)div  = 

hTCiZ 

Now  define 

(3.49) 


(3.50) 


«<■)  - 


(3.51) 


OLM  should  be  analytic  in  the  complex  plane  cut  £rom  0  to  oo  (the 
"cut  plane"),  vanishing  as  |z  |  -»  oo.  (3.48)  becomes 


aV;/\»-«»/A7u)=0  u>o 


(3.52) 


To  solve  this,  we  construct  a  function  X  (z)  which  is  analytic,  and, 
along  with  its  boundary  values  X^(u),  non-vanishing  in  the  cut  plane, 
such  that 

X+(u)  _  .A^ 


(3.53) 


Consider  the  function 
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(3.54) 

In  (  A^/  A  )  i*  that  branch  which  vanishes  as  u  — ^oo.  Then 
^  (z) 

e  satisfies  (3.53),  since 

(«'’)  _  .4.I.VA-]  A* 

[Tr  ■  "A- 

(3.55) 

However,  since  In  [  A  ^(0)/  a”(0)]  *2iri,  we  find 


(3.56) 


in  the  neighborhood  of  z  =  0.  Thus  an  appropriate  X(z)  is  given  by 


X(i) 


(3.57) 


Combining  eqs.  (3.52)  and  (3.53),  we  find 


XVjqVv)  = 


(3.58) 


which  shows  that  X  (z)  a  (z)  must  be  an  entire  function.  But  since 


a(K)  =  o 


(3.59) 


■  64. 


X(z)  ^  (z)  must  be  a  bounded  function.  Thus 


0.(2)  = 

AC2/  (Liouville's  theorem)  (3.60) 

We  may  obtain  the  index  of  refraction  from  this  general  result 
s  imply  by  noting  that 


r^v^A(v)dv  ^  . 

!  f  **  P^A(v)clv  (  I  \ 


X'U) 

X^2> 


Thus 


71  s=. 


cjfvAlv)olv 
J  P*A(v)alv 

o 


-CXVq) 

X(o) 


This  may  be  simplified  by  means  of  the  identity 
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(3.61) 


(3.62) 


X(z)X(-H)=^Afe) 

(3.63) 


which  is  most  easily  verified  by  noting  that  in  view  of  (3.  53),  R(z) 
defined  by 


A(2) 

-  (.  XWX(-2)  (S-SO 
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is  an  entire  function.  But 


^(2) 

!»»-*•  w 

Hence, 


/2i  ->«  -6  a*- 


(3.65) 


R(2)-I 


(Liouville 'a  theorem)  (3.66) 


It  follows  that 


X(o)  = 


(the  phase  is  easy  to  verify).  Thus 


(3.67) 


n  *  x*(o)/6 


(3.68) 


which  is  jvist  the  zero  temperature  result  with  the  correction  factor 


X'  (0).  X'  (0)  may  be  approximated  through  the  use  of  the  identity 

.00 
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Xfz)  *  X(o)  f  2  t 


_r/\v 

u(u- 


ru)  -A7t*) 


2)X(-U) 


clu 


O  (3.69) 

To  prove  this  second  identity  we  make  use  of  the  fact  that  [X(z)  -  z  -  X(0)  ]  z*** 
is  analytic  in  the  cut  plane,  vanishing  as  |  z  |  oo.  Hence,  by  Cauchy's 
theorem: 


_L  f  x^u)-x“^u) . 

X  AtiJzVa'-z)  14  (u -2) 


:(*)-? -xro; 


(3.70) 
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But 


(3.71) 


from  which  the  identity  follows. 


Using  (3.69)f  we  obtain  the  result 


X'(o)  =  I  + 


I 

ixii 


i 


1*“^  X(-^) 


du 


/  -  -L  fdt  f^'uF(u)clu. 
€  U>*J  X(-^) 

9 


(3.72) 

-mu*  /2  d 

Since  F(u)'ve  .  it  is  reasonable  to  approxinnate  X(-u)  »  X(0) 

in  the  integrand.  Thus 


XVo;  /  - 


j _ ^ 

ZeKM 


1 1 


«()*  Ifmc* 


(3.73) 

to  first  order  in  (2G/mc*)^^^.  Higher  order  corrections  tnay  be 
obtained  by  iteration. 
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7.  Digression  on  the  Outgoing- Wave  Boundary  Condition 

In  treatnnent  of  half-cpace  problems  in  this  chapter,  we  have 
somewhat  arbitrarily  reqtiired  that  there  shall  be  only  "outgoing  waves" 
(section  4)  or  that  the  frequency  should  be  considered  to  contain  a  small 
positive  imaginary  part  (section  5).  The  second  prescription  actually 
incorporates  the  first,  and  both  are  the  results  of  attempts  to  correct 
the  mistake  made  in  posing  the  problem  of  a  wave  which  has  been  in 
existence  for  an  infinite  length  of  time,  rather  than  treating  the  'actual 
build-up  of  the  wave  as  an  ixutial  value  problem. 

By  way  of  illustration  we  solve  the  initial  value  problem  for  the 
zero  temperature  plasma,  a  case  where  the  results  are  well-known. 

The  plasma  fills  the  half-space  z  >  0  and  for  t  <  0  the  electromag¬ 
netic  fields  are  confined  to  a  fixute  region  on  the  vacuum  side.  In  general, 
for  x-polarized  disturbances 


H<o 


(3.74) 


Taking  the  Laplace  transform  with  respect  to  time,  we  obtain 


(3.75) 
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We  need,  o£  course,  the  constitutive  relation  between  D  and  E.  As 
a  rule  one  writes 


(3.76) 


The  rationale  for  this  is  well  known.  If  we  consider  a  disturbance 

Ci>  t 

of  time  dependence  e  ,  then  the  equation  of  motion  for  each  electron. 


^  i  -a  g  £ 

r  =  -  ^ 

m 


(3.77) 


automatically  yields  as  a  solution  for  the  polarization 


p= 


(3.78) 


and  hence 


"*  **  -a  /  COm  \ 

V  =  E-^4icp 

(3.79). 

This,  of  course,  assumes  that  the  oscillation  was  initiated  far  enough 
in  the  past  for  the  polarization  to  build  up  to  its  final  value. 

In  the  initial  value  approach  we  write 


(Euler's  identity)  (3.80) 
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5  a)  =  Eft)  +«'ne?(t) 


=  E(f)+ 


and 


(3.  81) 


(3.  82) 


where  we  have  used  the  convolution  formula  to  write 


|e'/’‘’dt[ft-f)Eft')dt'  = 

0  (3. 83) 

Thus  the  expression  (3.  76)  is  valid  if  it  is  interpreted  as  the 
complex  Fourier  time  transform  of  £q.  (3.  81)  where  we  have  the 
usual  correspondence  u  =  ip.  As  an  indication  of  what  we  can  expect, 
we  note  that  for  p  on  a  typical  Laplace  inversion  contour,  (parallel  to 
the  imaginary  axis  in  the  right  hand  plane)  u  will  have  a  positive  imag¬ 
inary  part. 

In  order  to  avoid  loss  of  generally,  we  shsdl  write  instead 
of  (3.  82) 

Dp  =  tpEp 

(3.  84) 

We  now  take  the  Fourier  transform  of  (3.  75)  with  respect  to  the 
spatial  variable  z. 
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-(kE,(k)  =  -i{pBfW-B(k,o)} 

-  iU  Bp(k)  =  -  -i(k,o)} 


(3.  85) 


where  the  Fourier  transform  is  defined  by 


Ep(k)  = 


(3.  86) 


Eliminating  B^(k),  we  obtain 


+  ^Jp(k)  =  B  E(l<,»)+  tBCK.o)  =  i(k-!t)E((c^o) 


(3.  87) 


where  we  have  used 


B(k,(i)  =  E(k,0) 

( 

which  represents  the  fact  that  the  initial  disturbance  consists  of 
waves  traveling  towards  the  interface  with  velocity  c.  Now  make 
the  separation 

ip(k)  -  Ep(k)  +ip(k)  =  +/jEpf»je‘''*dz 

D,(k)  +Ep(k) 

^  ( 

This  yields 

(k*4  ^)Ep(k)  +  (k'4  ^‘)E;(k)  =±(k.i^)g^k,o) 


(3.  88) 


(3.  89) 


(3.90) 


I 
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Be sides  enabling  us  to  express  O  conveniently  in  terms  of  E, 
the  separation  (3.  89)  places  our  equation  in  a  form  suitable  for  treat- 
ment  by  standard  complex  variable  techniques.  The  function  Ep(k) 
may  be  analytically  continued  into  the  upper  half  k-plane,  and  is,  in 
fact,  analytic  throughout  the  upper  half  plane  and  vanishes  as  k-^i  a>. 
Similarly,  (k)  is  the  boundary  value,  as  k  approaches  the  real 
axis  from  below,  of  a  ftmction  analytic  throughout  the  lower  half  k-plane, 
vanishing  as  k — i  oo.  Since  the  initial  disturbance  is  confined  to  the 
region  z  <  0,  £  (k,  t  s  0)  is  a  boundary  value  in  the  same  sense  as  £^  (k), 
and  the  (  -  )  superscript  on  the  right  hand  side  of  £q.  (3.  90)  serves  to 
remind  us  of  that  fact. 

Equation  (3.90)  is  of  the  form 

which  may  be  solved  by  the  same  technique  used  in  section  8.  We  must 
factor  the  function 


A^k) 

X*(K) 

-X'(k) 

(3.  92) 

where  (k)  is  the  boundary  value  of  a  function  analytic  and  nonvanishing 
upper 

in  lower  plane.  In  this  case,  the  resxxlt  is  obtained  by  inspection: 

(3.  93) 
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(It  is  assumed  that  we  may  define  the  square  root  in  such  a  manner 
1/2 

that  (p  c  lio>  io  the  right  half  plane  whenever  p  lies  in  some 

appropriately  chosen  half  plane  Re  p  >P^> )  Then 

xN  k)E;  (w)  -  w  =' &(t)  =  i 


from  which 


>  fl  (■.'-‘•f)  i-(lc',o)dk' 


(3.  94) 


95) 


or 


(3.  96) 


Inversion  of  the  Fourier-LaplKce  transform  is  straight-forward: 

Ufcl_Trv^,E(-«o)e  e  du, 


(z-^o) 


(3.97) 


'  (H.O) 


where 


f(U3)=  S_; 


liO 


(the  dielectric  coefficient) 


(3.  98) 
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and  where  we  have  written  throughout  to  indicate  that  whenever 
u  appears  in  the  integrand,  it  should  be  thought  of  as  having  a  small 
imaginary  part  which  is  allowed  to  tend  to  zero  after  the  integration 
is  performed. 

We  interpret  the  different  terms  in  E  (z,  t)  as  superpositions 
of  transmitted,  incident  and  reflected  waves  respectively.  Having  done 
this  we  see  that  the  transmitted  and  reflected  electric  field  amplitudes 
for  waves  of  a  given  frequency  are  given  in  terms  of  the  amplitude  of 
the  incident  field,  viz; 


* 


— - £.  =  -i:—  E 

I  -t-  *•  l-tnfu))  ‘ 


p  ^  F  -  /  -  p, 

^  “  T-t  ‘  ~  ‘ 


(3.99) 


(n{b>)  is  the  index  of  refraction) 

Note  that  the  solution  to  the  initial  value  problem  automatically  satisfies 

the  condition  that  the  electric  field  inside  the  plasma  does  not  contain 

incoming  waves.  This  condition  is  usually  introduced  ad  hoc. 

We  have  shown  that  the  initial  value  problem  has  a  unique  solution. 

without  the  imposition  of  arbitrary  subsidiary  conditions.  In  fact,  we  have 

constructed  the  solution.  This,  therefore,  is  the  "correct"  way  to  solve 

the  problem,  or  rather,  it  is  the  correct  way  to  pose  the  problem.  What 

then  is  meant  by  a  solution  to  the  artificial  problem,  which  specifies 

*  i  u  t 

steady-state  time  dependence  e  ?  Nothing  more  than  the  interpre¬ 


tation  of  the  above  result  as  a  superposition  of  steady-state  solutions. 
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Working  backwards,  we  will  agree  that  we  have  a  correct  solution  for 
the  reflected  and  transmitted  amplitudes  E^(u),  E^(u).  if  a  superposition 
of  such  solutions  yields  the  correct  solution  to  the  well-posed  initial  value 
problem.  This,  of  course,  is  a  general  criterion  which  should  be  applic¬ 
able  to  the  case  where  the  medium  fills  a  slab  rather  than  a  half- space, 
or  even  if  the  mediiun  does  not  satisfy  phenomenological  equations  of  the 
form  (3.74). 


The  above  result  (3.  97)  suggests  the  assertion  that  in  considering 
more  general  problems  from  the  steady-state  (time  dependence  e~  ^  ^  *) 
point  of  view,  a^  ambiguities  will  be  correctly  resolved  by  assumine  that 
u  has  a  small  positive  imaginary  part  and  requiring  that  the  fields  vanish 
as  z—*CD,  Of  ^purse,  in  any  particular  practical  case,  one  would  use 
this  prescription  to  solve  for  E^(u)  and  E^(u),  and  then  justify  it  a 
forteriori  by  examining  the  solutions  so  obtained  in  the  light  of  the  corre¬ 
sponding  initial  value  problem. 

The  preceeding  remarks  may  be  seen  in  a  more  general  light  by 
understanding  that  causality  lies  at  the  root  of  all  these  considerations. 
Suppose  we  restrict  pur  attention  to  the  vacuum  side  of  the  interface,  with* 
out  any  restriction  on  what  lies  beyond,  except  that  it  be  a  linear  system. 
Then,  quite  generally 


(3.100) 
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Mathematically,  *«*’>*'  represents  a  disturbance  which  is 

riffKt 

traveling  to  the  [  ]  without  change  in  spatial  form.  The  condition 

that  the  initial  disturbance  consists  entirely  of  waves  propagation  towards 
the  interface  from  the  left  and  is  localised  to  the  left  of  the  interface 
requires  that  at  t  =  0  the  function  Fj  (s  -  c  t)  is  non-vanishing  only 
for  z  <0,  and  the  function  F2  (z  +  c  t)  is  non- vanishing  only  for  s  >  0. 

(Of  course,  the  expression  (3.  10)  is  valid  only  for  z  <0. )  Hence,  E.(u) 
and  E^(u)  are  boundary  values  of  functions  analytic  in  the  upper  half 
plane. 

As  before,  we  shall  interpret  the  functions  Fi  and  F2  as 
incident  and  reflected  fields,  the  reflected  fields  being  caused  by  the 
response  of  the  medium  behind  the  interface  to  the  incident  field.  Caus¬ 
ality  then,  is  simply  the  requirement  that  there  can  be  no  reflected  field 
until  the  incident  field  distrubance  has  reached  the  interface.  Moreover, 
in  this  interpretation,  since  the  medium  is  linear,  the  reflected  field  F2 
must  depen  d  linearly  upon  the  incident  field  F2  which  causes  it.  Thus 

‘  (3.101 

where  ^(u)  the  "response  function"  depends  only  upon  the  properties  of 
the  medium,  and  is  the  boundary  value  of  a  function  analytic  in  the  upper 
half  plane. 

All  this  serves  as  a  justification  for  including  a  small  imaginary 
part  in  u  when  solving  the  steady- state  problem- -time  dependence  ** 
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Since  all  field  ai^plitudes  are  boundary  values  of  analytic  functions, 
we  may  seek  steady- state  solutions  with  u  in  the  upper  half  plane, 
with  the  assurance  that  the  resulting^(u)  will  be  the  analytic  contin¬ 
uation  of  the  desired  response  function.  In  other  words,  we  insure 
that  the  field  amplitudes  will  possess  the  required  analytically  properties 
by  working  in  a  region  of  the  u  plane  in  which  they  are  known  to  be 
analytic. 


8.  Propagation  Along  a  Constant  Magnetic  Field 

If  the  plasma  is  located  in  a  steady  magnetic  field  (not  assumed 

"small*|,  then  the  Vlasov  equation  takes  on  an  added  degree  of  complexity- - 

in  particular,  it  becomes  a  differential  equation  in  the  velocity  variable. 

The  general  treatment  may  be  approached  through  the  "Fessel" 

19 

transformation,  which  transforms  away  the  magnetic  field  term  at  the 
start.  However,  in  the  special  case  of  wave  propagation  along  the  steady- 
state  magnetic  field,  the  problem  becomes  almost  identical  to  the  case  of 
zero  field.  We  sketch  the  results  here. 

In  the  linearization  of  (1.  34)  we  require 


"small" 


E 


"small" 


(3.  102) 


as  before.  On  the  other  hand,  now; 


3 


B(rA) 


"small" 


(3.  103) 
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1 


r 

i. 


whore  k  is  a  unit  vector  in  the  z-direction.  Note  that  if  F  (u) 
is  a  function  of  energy  alone,  it  will  be  independent  of  since 
there  is  no  magnetic  potential  energy.  Hence,  if  is  isotropic 
in  the  zero  field  case,  it  will  also  be  isotropic  in  this  case.  The 
linearized  equations  are  then 

+  tl-vf  t 

Tit 

Vx£  if  A)  -  -  i  <W1  ^ 

V,6  =  ^  V-6  -  o 

We  have  used  Eq.  (1.41)  to  simplify.  Otherwise,  except  for  the 
term,  these  equations  are  identical  with  Eqs.  (1.39).  Because 
of  the  presence  of  the  term,  the  Vlasov  equation  does  not  separate 
into  plane  polarized  modes.  The  equations  corresponding  to  (1.47x) 
and  (1 . 47  y)  are 


where  we  have  introduced  £7  =  -  eB^/mc  (the  cyclotron  frequency). 
However,  these  do  separate  into  right  and  left  circularly  polarized 
modes:  writing 
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]■ 


r 

I. 


E* 

B. 

(3. 106) 

9±  =  5x  i  98 

we  get 

/ 

.ul  3.  -i^9*  = 

^  I 

right 

circularly 

(3.107+) 

polarized 

modes 

-  p  -  i.  Is 

^  c'it  ^  t  J 

left 

circularly 

polarized 

(}t3-  * 

(3.107-) 

modes 

1  1  c  -  t  'JB. 

la  -  "  c 

-'  li  "  c  ^f^-*^** 

The  plasma  oscillations  are  the  same  as  in  the  zero  field  case,  since 

electrons  undergoing  longitudinal  vibrations  (  ||b  )  do  not  experience 

o 

any  "u  x  force. 

The  normal  modes  may  be  exhibited  by  reproducing  the  zero  field 
treatment  of  chapter  ZI,  in  fact,  the  remaining  equations  in  this  section 
will  be  seen  to  be  almost  identical  with  corresponding  equations  in  chapter  II. 
Fourier  analysis  in  the  time  variable  yields  (the  z  >  dependence  and 


the  ±  subscript  is  suppressed) 
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]' 

i. 


o 
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/ 

0  1  o 

2. 

E 

j 

?2 

\o  0  1  J 

w 

\ 

c  J 


o 


T- 
*  c. 


(3.  108) 


or: 


The  corresponding  adjoint  equation  is: 


O  o\ 
o  1  0 


o  o 


gV)  eJ  2^ 


(3. 109) 


I 


0 

0 

Tf 

*  t 


1(3.  110) 


or: 


•^2 


The  solutions  of  (3. 108)  and  (3. 110)  with  space  behavior  of 


(3.111) 


i  u/v  z 

the  form  e  will  satisfy  the  orthogonality  relation  (2.  15).  In 

this  case,  discrete  solutions  will  exist  whenever  v,  is  a  root  of 

J 

AM. 


-eo 


(3. 112) 


where 


.  n 

*  UJ 


(3,113) 
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The  continuum  ii^cludes  the  entire  real  axis.  The  solutions  to  (3. 


are 


X  4ttCI>C 
CO 


^xiev^ 

CO 


3  FM 

u-ecv 


for  the  discrete  modes.-  For  the  continuum  modes 


(U)  =  V*  P  +Xl’J)Slu-cev) 

X(u)  =  ±  [a’'(w) 


The  adjoint  solutions  are 


x  y)ei/^c 
moj 


pt  _  'Tfievc’’ 

V  mitu 


i; 


3j _ 

U  -0{v 


for  the  discrete  modes,  with 


X(v)F(MV)  =  Mv) 


for  the  continium  modes.  By  direct  calculation,  we  find 


N(v)=^|V 


X‘(v) 

F(«v) 


108) 


(3.114) 


(3.  115) 


(3,  116) 


(3.117) 


(3.  118) 


which  may  be  written  in  the  form  (cf.  (2.  85)) 
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-L  1  ^ 

N(v)  2Ki 


J _ 1 

/\\v) "  A"(v) 


} 


(3.  119) 

The  discrete  spectrum  also  is  modified  by  the  presence  of  the  magnetic 
field.  Consider  first  the  or  "ordinary"  modes, 


(3.  120) 


At  zero  temperature  (  F  (u)  »  5  (u)  )  we  find 


v‘ 


I  - 


o): 


(jdCoJ'Irfl) 


1. 

=  C) 


which  yields 


(3,  121) 


a)  nq  complex  roots:  two  real  roots,  when  w^(Wp  +  (/i/2)*)^^^  -/1/2 

1/2 

b)  two  complex  roots  at  v=  ±  i  ^  when  w<  (w^  +  (JX/Z)*)  0,/Z 


with 


.-i 


_,v 

At  finite  temperaturss,  the  situation  is  similar.  There  will  be  |  ^/discrete 
,  1/2 

modes  when  b)^(u  *  *  )  -  Cl/Z.  The  real  roots  have  wandered  off 

the  real  axis  to  the  other  sheet  of  the  double- valued  function  ^(v)-  Zn  other 
words,  if  we  consider  the  expression  (3,  120)  for  A(l^)  with  v  in  the  upper 
half  plans,  we  will  find  no  zeroes,  but  if  we  analytically  continue  it  into  the 
lower  half  plane,  we  find  (among  othe.'s)  two  zeroes  which,  at  low  temperatures 
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at  any  rate,  are  near  the  reel  axis.  Of  course,  these  seroes  do  not 
represent  actual  discrete  modes  because  the  expression  (3.  120)  for  A  (  v) 
does  not  represent  its  analytic  continuation  across  the  real  axis.  This  is 
no  different  from  the  case  when  there  is  no  magnetic  field. 

The  or  "extraordinary"  modes  are  somewhat  different.  At 
zero  temperature 

^‘(l  = 

(3. 123) 


I.e. 


a)  no  complex  roots:  two  real  roots,  when  u»>(Up*  +  Cft/2)*)^^^  +il/2 

b)  two  complex  roots  at  v  r  i  i  v.  _  when  SI  <  w  <(u*  +  Ul/2)*)^^^  +S1/Z 

P 

with 


(3.  124) 


c)  two  real  roots  when  u>  < 

In  cases  a)  and  b),  the  situation  is  the  same  as  for  the  ordinary  mode.  There 
will  be  0  and  2  discrete  modes  for  cases  a)  and  b)  respectively.  In  case  c), 
something  new  develops.  The  real  zeroes  of  the  zero-temperature  A(  v) 
become  complex  at  finite  temperatures  and  give  rise  to  four  complex  roots, 
two  each  in  the  upper  and  lower  half  planes,  of  the  dispersion  relation 
A  {  v)  -  0  with  corresponding  discrete  modes.  In  fact,  if  we  were  to  apply 
the  argument  principle  to  obtain  the  number  of  roots  in  (say)  the  upper  half 
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plane,  we  would  find  that  the  curve  corresponding  to  figure  2,  p.  44,  would 
look  like 


which  encircles  the  origin  twice,  so  A  argA  =  4  ir .  The  discrete  roots 
are  symmetrically  located  at  ±  v^,  and  at  low  temperatures,  near 

the  zero- temperature  values 

-'4 


^  I  I  + 


co; 


co((A)+X)))  C 


(3.125) 


We  find 


with 


+  T  lri«H 


eo 


'  ^  w*  I  +. 


•4 


(3. 126) 


(3. 127) 


(3.  128) 
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The  presence  of  these  discrete  modes  has  led  to  some  confusion  in 
the  literature.  For  example,  in  the  transmission  of  electromagnetic  rad¬ 
iation  through  a  half-space  (in  the  extraordinary  mode  with  CO<Jl ),  the  spatial 
dependence  of  the  amplitude  far  from  the  boundary  will  be  given  by  the  slightly 
damped  discrete  modes.  In  the  paper  cited  above,  this  behavior  is  incorrectly 
attributed  to  Landau  damping,  a  phenomenon  discussed  in  the  next  chapter. 


Chapter  IV.  OTHER  APPLICATIONS 


1.  Initial  Value  Problema 

In  thf  initial  value  problem,  the  boundary  value  a  specified  are  the 
values  of  all  field  quantities  for  all  points  in  space  at  a  fixed  initial  time,  t^. 

It  is  desired  to  deternoine  the  complete  behayior  of  the  system  for  later  times. 
There  is,  then,  a  basic  difference  between  the  initial  value  problem  and 
boundary  value  problems  of  the  kind  considered  previously,  where  the  data 
was  specified  for  all  times  at  a  finite  number  of  space  points.  The  nature  of 
the  initial  value  problem  boundary  conditions  requires  that  the  problem  be 
Fourier  analysed  in  the  space  variable  rather  than  in  the  time  variable  as 
previously  (cf  eq.  (2. 1)).  Also,  as  a  rule,  in  the  initial  value  problem,  the 
plasma  is  taken  to  be  an  infinite  medium.  Nevertheless,  in  terms  of  the 
resulting  equations  to  be  solved,  these  differences  are  purely  formal.  Actually, 
the  singular  eigenfunction  method  we  have  been  using  has  been  applied  more 
often  to  the  initial  value  problem  than  to  other  types  of  boundary  value  prob¬ 
lems.  In  fact,  the  initial  value  problem  was  the  subject  of  van  Kampen's 
paper^  in  which  singular  normal  modes  were  introduced. 

2.  The  Initial  Value  for  Transverse  Modes 

This  is  by  far  the  most  often  treated  problem  in  plasma  oscillations. 

The  discussion  in  this  section  is  included  as  background  material,  and  contains 

ip  — —  yj 

Montgomery  and  Gorman  have  shown  that  in  certain  cases,  slab  geometry 
may  be  treated  by  using  the  method  of  images  to  convert  the  slab  problem  into 
an  infinite  medium  problem. 
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only  an  outline  of  the  reaulti  in  the  quoted  papers. 

Consider  the  spatial  Fourier  decomposition  of  eq.  (l.dTs),  I.e. 
let 

*  (4.  1) 

Then 

i9(V)  + 

OX  (4. 2) 

and 


ikE 


4irp  =  4ir€f  gdu 


(4.3) 


The  coupled  equations  (4.2,  3)  were  first  solved  by  Landau^^  using  the  Lqplace 
transform.  The  result  for  the  Laplace  transform  of  the  electric  field  is 

4' 

21 

Van  Kampen  utilised  the  normal  mode  approach.  Under  the  ansatz  that 
the  solutions  have  time  dependence  e’  ^  ^  eqa.  (4.2|3)  reduce  to 


f  rwdt  = 


^  r" 


Cu-v)g;;(u)  -  qju.)cjiu 

"  T^l  9t,(vi)du 


(4.5) 
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The  normal  modea  are  then 


+  X(v)S(u-v) 


(all  real  V ) 


with 


X(v)  =  I  ( ^Uv)  +  Alv)]  =1-^1  Pf-^^ 

"to 


(4.6) 


(4.7) 


Discrete  modea  would  also  exist  for  those  v  which  are  roots  of 


A(v)  =\-^. 


''""FW)dM 


=  O 


-60 


(4.8) 


but  for  the  case  when  (Ti)  is  Maxwellian,  the  discrete  spectrum  is  empty. 
Note  that  the  denominator  in  eq.  (4.4)  is  just  A  (’^). 

Van  Kampen  proved  that  the  g^(u)  formed  a  complete  set,  by  means 
of  the  singular  integral  equation  technique  discussed  in  chapter  U.  He  proved 
completeness  by  constructing  the  solution  to  the  integral  equation  implied  by 

oo 

Uv;cj^(u)clu 

(4.9) 

where  g  (u,  o)  is  an  arbitrary  function,  (the  initial  perturbation  in  this 
instance).  The  solution  is 


pOO 

g(u^o)=s.  A(v;g^(u)clu 
*>  J-oo 


A  -  X  -  -  6r'(v) 

A"(v) 


(4.  10) 


where 


Ziti  Lpo  ti  -z 


(4.11) 


r 
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j 

The  tolution  to  the  initial  value  problem  ia  then  aimply 

1 

t 

£(tj=  Arv)e-‘'“'''dv 

*'-00 

(4.  12) 

i 

14 

Case  placed  the  normal  mode  treatment  on  a  more  formal  basis. 

1 

He  showed  that  the  g  (u)  are  orthogonal  to  the  solutions  of  the 

V 

adjoint  equation: 

(u-v)a^(u)=  ^  F'(w)0t(u)c|u 

■'  J-oo 

(4.  13) 

That  is 

1 

( 

i 

?* 

.00 

1  *  U(v)8iv-v') 

J_00 

(4.  14) 

1 

The  adjoint  solutions  are 

7- 

1 

1  • 

^  +  A(v)S(«-v; 

(4.  15) 

with 

F'(v) 

(4.16) 

Then  an  arbitrary  initial  perturbation  g  (u,o)  maybe  expanded  in  terms  of 

the  (u)  with  the  expansion  coefficient  given  by 

I  - 

r  ^ 

A(v)  =  j;^)J  9>)9(u,o)dM 

••  Oo 

(4.  17) 

! 

! 

The  normalization  coefficient  is  found  to  be 

ji' 

li: 

I  T 

N(v)  =  ^  ^ 

(4. 18) 

1  I 
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which  may  be  written 


I 


1 


N(v)  Aici  % 

kr 


JL  l-_ 

Alv)  Mv) 


(4.  19) 


Case  also  pointed  out  that  in  view  of  the  Plemelj  formulae,  these 
apparently  different  approaches  to  the  initial  value  problem  are  actually 
equivalent.  That  is,  it  follows  almost  by  inspection,  that  the  normal  mode 
expansion  (4.  12)  for  E  (t)  t  >0  and  the  Laplace  inverse  of  (4.4)  are  actually 
the  same  contour  integral  taken  along  two  equivalent  paths. 

This  last  result  is  quite  general,  and  is  valid  even  when  the  steady- 


state  distribution  function  ^  (u)  gives  rise  to  discrete  modes. 

o 


3.  Stability  and  Landau  Damping 

It  is  clear  that  if  there  is  a  complex  root  of  the  secular  equation 

(4.  8),  the  system  will  be  unstable,  since  (at  least  for  some  initial  perturba- 

•  i  k  V*  t 

tions)  the  electric  field  will  contain  a  term  e  o  which  grows 

exponentially  either  for  k  or  -  k.  It  has  been  pointed  out  that  for  Max¬ 
wellian  (u)  there  are  no  discrete  roots.  In  fact,  it  is  siinple  to  show  that 
the  discrete  spectrum  will  be  empty  if  F  (u)  represents  any  "single-humped" 
distribution.  I.e. 


p'(u)  ^  O  U  ^  “Mo 


(4.20) 
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7’ 

I. 


We  prove  this  by  applying  the  argument  principle  to  A  (z) 
traverses  the  contour  C 


The  contour  C 

I  Along  the  large  semicircle 

?  I* 

j  1.  AU)  ^  I 

[ 

I 


as  z 


(4.21) 
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Along  the  real  axia 

A*(x)  =  I  -  -Tri^>'(x) 

and  since  F  (u)  is  single-humped 


-40 


(4.22) 


Jry)  A^(^)  ^  0 

X  ^  14^ 

(4.23) 

fs  4 

Thus  the  image  of  the  path 

L4--»iJ 

lies  entirely 

below 

Labove, 

the  real  axis. 

The  image  of  the  path  3  -  4  •  1  will  then  cross  the 

real  axis  only 

once,  at 

i 

the  point  A 

i‘ 

i 

7* 

Km  =  I  ~  ^ 

f  FV“)d»  >1 

)  U-Uo 

(4.24) 

1. 

since  is  negative  definite. 

Thus  the  image  of  the  contour 

C  looks 

u  -  u. 


something  like 


which  does  not  encircle  the  origin,  and  there  are  no  discrete  roots. 
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When  there  are  no  discrete  modes,  the  system  will  be  stable.  Thia 
follows  from  the  normal  mode  representation  (4.  12)  of  the  electric  field. 

For  any  g  (u,  o),  the  function  G  (z)  will  be  analytic  in  the  cut  plane,  and 
if  (as  we  will  assume)  g  (u,  o)  is  reasonably  well  behaved,  the  boundary  valubs 
(V  )  will  be  non-singular.  Thus,  for  t  >  0 


E(t) 


- 


A^fv) 


(4.25) 


(the  term  involving  — gives  no  contribution,  since  —7  is  the  boundary 

A  (  v) 

value  of  a  function  analytic  throughout  the  lower  half  plane). 

The  asymptotic  behavior  (t— »oo)  of  E  (t)  is  governed  by  the  singularity 

G+ 

nearest  the  real  axis  of  the  analytic  continuation  of  — —  into  the  lower  half 
plane.  Landau  showed  that  if  ^  (u)  is  Maxwellian,  the  integral  expression 
(4.  8)  with  Pmz  >  0  when  continued  into  the  lower  half  plane,  has  a  zero  (for 
small  k)  given  by 


CO^  =  0)^^  +•  I  Uloi 


(4.26) 


where 
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Jackson  has  shown  that  a  slightly  better  approximation  is  given  by 

U),r  «  +3k‘^ 


(4.28) 


At  any  rate,  a  non-singular  initial  disturbance  gives  rise  to  an  electric 
field  which  after  long  times  undergoes  damped  oscillations or  ^  e  ‘^ol^. 
This  "Landau"  damping  may  (in  keeping  with  the  van  Kampen  singular  normal 
mode  picture)  be  thought  of  as  the  result  of  the  phase-mixing  of  an  initial 
disturbance  composed  of  oscillations  which  propagate  with  a  continuous  spectrum 
of  velocities. 


4.  Landau  Damping  of  Transverse  Modes 

The  initial  value  problem  for  transverse  modes  may  be  treated  in  the 
same  manner.  The  discrete  mcdes  of  the  system  would  have  space-time 
behavior  e  where  the  are  the  zeroes  of  the  dispersion  func- 

tion  ^  (  v).  In  this  case,  however. 


reo 


'Leo 


(4.29) 


with  fixed,  real  k  and  complex  v  ^  w/k,  and  it  is  easy  to  show  that  there  are 
no  discrete  modes.  Hence,  for  positive  ^negative]]  k  the  long-time  behavior 
of  the  system  is  of  the  form~e^^*  where  is  the  zero  lying 

nearest  the  real  axis  of  the  analytic  continuation  of  A  ^( '')  j^A  ”  (  '^)~|  into 
the  lower  ^pperj  half  plane. 
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We  note  that  the  analytic  continuation  of  (4.29)  from  the  upper  half 
plane  into  the  lower  half  plane  is  given  by 


Ac„t  -  ^*«ivRv) 


(4.30) 


Then,  if  we  assume  that  k  lies  near  to  its  zero-temperature  value  of 
2  221/2 

(co  +  k  c  )  ,  we  may  compute  k  v  by  means  of  successive  approx- 

P  ° 

imations.  In  particular,  the  attenuation  decrement  Im  (k  ia  found  to  be 

e'-'" 


31) 


5.  Longitudinal  Oscillations:  Fixed  Frequency  Modes 

The  longitudinal  oscillations,  like  the  tranverse  oscillations,  may  be 
analyzed  into  fixed  frequency  modes: 

•  J 

iu>t 


E/i,t)  =E(a)e 

The  equations  corresponding  to  (4.2-3)  are 


(4.  32) 


or 


lU) 


\ 


E(i) 


nlr 


=  Hf 


The  adjoint  equation  is 


a  YV  =tV 


n 


(4.33) 


(4.  34) 


(4.35) 
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where 


(4.36) 


Here  we  must  modify  the  procedure  used  in  the  transverse  mode 
case.  Zero  is  an  eigenvalue  of  the  operator  H.  I.  e. 


HT.  =  0  =  y;h+ 


(4.37) 


with 


_  /  Wiu)  “  ' 


1 


(4.38) 

Thus  the  set  of  solutions  of  (4.  34)  contains  the  spatially  uniform  solution  ^ 

•  i  it)  S  /  V 

as  well  as  the  modes  obtained  by  the  separation  ^  e  These  last 

solutions  are  pbtained  in  exactly  the  same  manner  as  the  transverse  modes 
derived  in  chapter  U. 

The  continuum  modes  are 


'lip 


(a  |ix)  =  ^  ■hA(v)^(u-v) 

’  60^  ^  U-V 


i  00 


(4.39) 


t 


p— :  ■+■  A(v)^ru-y) 

(Aj^  Ly 


oo 

_  naiv 


(4.40) 
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with 


-Xiv)  =  X(v)F'(v)=  I  - 

x(-v 

-00 


(4.41) 


Whenu  ^  u_,  there  are  no  roots  of 
P 


A(i)=  I  - 

'^-00 


(4.42) 


and  when  u  <  u^,  there  are  two  imaginary  roots  1  i  and  corresponding 
discrete  modes.  Note  that 


A(»)  =  I  - 


Ab  before,  we  have  the  orthogonality  relation 


.i,  ,  ,,[ir'^VF'iv4i  +  A‘'(W 

N  v)  =v-— - 

F(v) 


which  may  be  written 


where 


(4.43) 


(4.44) 


(4.45) 


N(v)  AV)  A"(W 


(4.46) 


Of  course,  the  discrete  modes,  if  any,  are  orthogonal  to  the  continuum  modes 


and  to  each  other. 
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The  set  (discrete  and  continuum)  is  complete  when  supplemented 

by  the  spatially  independent  solution  (J;  The  completeness  proof  is  straight¬ 
forward.  We  attempt  to  expand  a  state  function  in  terms  of  the  continuum 
modes . 


or 


oo 

=  I 

J-oa 

g(«)  -  F'M  P|  +AmA(u) 

4TCe  r‘ 


I  vA(v)olv 

'M 


Combining: 


(4.47) 

(4.48) 

(4.49) 


“U  (  0(u)  4-  •"  A'('*4)Cf(M) 

(4.50) 


where  we  have  defined 


n/2)  =  ±.  f"vA(v)dv 

(4.51) 

a(-)  is,  by  assumption,  analyww  in  the  cut  plane,  vanishing  at  least  as  fast 
as  l/|z|  as  |z|-.»oo.  In  fact. 


JU»\ 

Ixl  -♦•0 


-2(^(Z)=  Aim 

lai-^so 


X 

ZTCi  we 


(4.52) 


which  is  a  subsidiary  condition  that  ^(z)  must  satisfy  in  order  that  A  (  v) 
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given  by 


vMv)-  CC'i^v)  -  CC(v) 

shall  be  a  solution  of  (4.47). 

The  solution  of  (4.50)  is  clearly 

A(Z)i^  U-Z 


(4.53) 


(4.  54) 


which  is,  in  fact,  analytic  in  the  cut  plane  provided  that 


(cf.  eq.  (2.79))  (4.55) 


The  subsidiary  condition  (4.52)  beconrves 

.oo 

u 

t0‘ 

This  may  be  written,  using  (4.  38) 


too 

Hite 


-  — 5 ^  f  u  olu  (<3(u)  +  ^  EF  (**)) 


( 'vj .  ?  t)  =  0 


(4. 56) 


Thus,  any  function  which  is  orthogonal  to  ^  as  well  as  to  the 

1.0^  may  be  expanded  in  terms  of  the  continuum  I  ij;  1  •  This  completes 
t  i  L  v’ 

the  proof,  since  an  arbitrary  function  ^  may  always  be  written 


where  is  orthogonal  to  4*"^  and  to  the  4^^ 

o  ±i  V 


(4,57) 
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6.  Boundary  Value  Problems  (Longitudinal  OtcilUtion*) 

By  virtue  of  the  completenese  proof  of  the  last  section,  we  may 
use  the  representation 


-f  (j,u)  =  Tam  dv  +  Ajt 


(4.  58) 

and  apply  the  boundary  conditions  appropriate  to  the  problem  under  consid¬ 
eration  to  evaluate  the  expansion  coefficients.  We  have  assumed  CO  >  CO^ 
only  to  avoid  the  necessity  of  carrying  the  discrete  terms  along  in  the  calcu¬ 
lation. 

One  consideration  may  be  applied  to  the  distribution  function  at  the 
start.  Regardless  of  the  nature  of  the  reflecting  surface,  there  will  be  no 
net  current  carried  by  the  electrons  across  a  plasma  boundary.  Hence,  at 
a  boundary  z  s  a.  Jug  (a,  u)  du  must  vanish,  which  implies 


(-yf,  ?  f (i.A)  ^  j _ 

^  1-^ 

= 

I 


f  E(a) 

-CO 


II 

O 


On  the  other  hand,  if  the  plasma  fills  a  half-space,  then 

Ao  “Ef") 

which  follows  from 

EU)=  [  A(v)£yC'’'*clv  +A. 

J-  00 


(4.59) 


(4.60) 


(4.61) 
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since  the  integral  term  vaniahes  as  z  — ^a>.  by  virtue  of  the  Rienvann- 
L^esque  lemma.  Thus,  the  familiar  result 


EM* 


(4.62) 


for  a  half-space  with  boundary  at  z  *  0. 

Now  consider  the  boundary  condition  for  the  electron  distribution. 
If  the  boundary  is  perfectly  reflecting 


g(a,u)  =  g(a^-M) 

This  condition,  together  with  the  relations 

V'W)* 

(assuming  F  (u)  >  F  (-  u) ) 

F  «  —  F 
yields  ^ 


(4.63) 


(4.64) 


/  0  \ 

which  is  of  the  form  of  an  expansion  of  the  function  I  _  I  in  terms  of 

the  (I,  and  t|i  ,  and  the  expansion  coefficients  may  be  found  by  means  of 
V  o 

the  orthogonality  relations: 


2E(a) 


Nli^) 


(4. 65) 


Thus,  if  the  electrons  reflect  specularly  at  every  boundary,  we  have 
immediately: 

a)  For  the  half-space,  boundary  at  z  «  0 
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r 

i 


i 


I 

i 

I 


A(v)  -  Atw)  “  2  EM  -fig  jji 


A(v)-o  (v<o) 
A„  =  E(«)  = 


which  yields 


Eli)  e 


•(ut 


■  iiet 


{£ 


I 


Hlv) 

(outgoing  wavo  condition) 

(4.66) 


A(v)Eye‘''^dv+A,V 


(4.67) 


b)  For  the  slab 

A.  . 


E(o)  _  E(q) 
ITaT’  I  -  i4 


A(v)-A(-v)  -  A(v)c‘^‘'-A(-v)e“*®"=U(o)^iS 


(4.68) 


Thus. 


u) 

I  - 


I  _  ^  j55v; 


(4.69) 


and 


''  J  (4.70) 

Of  particular  interest  is  the  impedence  Z  of  this  "plasma  capacitor,  " 
defined  by: 


Z  j  -  feiDdi 


0 


(4.71) 
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r 

J. 


r 

1. 


where  J  is  the  net  current  flowing  into  the  condenser,  end  is  given  entirely 
by  the  displacement  current  at  z  «  0,  since  there  is  no  net  electron  transport 
across  the  boundary,  i.e. 


Hence 


J 

A 


^£(0) 


J 


_L_  i. 

iwCo  E(o)a  iu)C^ 


(4.72) 


(4.73) 


with 


C 

®  4ica 

the  capacitance  of  the  condenser  without  the  plasma. 


(4.74) 


^  ,  the  "effective  dielectric  coefficient"  may  be  written  in  terms 

^  eff 


of  the  function  A(z),  since 


[W=E(o,(t^+z5 


V  I  -  COST?  a)i 


J  (4, 


75) 


where  the  singularities  in  the  integrand  are  to  be  avoided  by  ascribing  to 
u  a  small  positive  part.  Thus,  we  may  apply  (4.46)  and  write  the  result 
as  a  contour  integral,  which  may  be  evaluated  in  a  residue  series: 


•eff 


(4.76) 


10*  odd  n 

Note  at  at  zero  temperature,  the  term  in  braces  vanishes  identically. 

Finally,  we  may  treat  the  half-space  with  diffuse  reflecting 
boundary.  In  this  case 

g(o,w  -  CE(o)F(«)  M>o 
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i 


which  represents  the  fact  that  the  wall  remits  the  electrons  in  a  thernnal 
distribution  at  the  same  temperature  as  the  plasma.  C  is  a  normalisation 
constant  to  be  determined.  (4.77)  and  (4.59).  along  with  the  outgoing  wave 
condition  A  (  v)  =  0  (v  <0),  imply 

Elo)  CrM  =  E(oi  |cF(u)  -  =| 

which  amounts  to  a  half-range  expansion  of  the  left  hand  side  in  terms  of 
half  of  the  continuum  modes.  This  is  of  the  same  form  as  the  expansion 
treated  in  Chapter  III.  sec.  6  .  and  we  only  sketch  the  results.  We  define 


vA(v)c|v 

2 


Then  (4.78)  combined  with  (4.59)  and  (4.61)  gives 


(4.79) 


E(o)  VI  [cF(m)  ~  V) A%)  -  C^"(m)A"(W 


U>0 


(4.  80) 

C\,(z)  is  to  be  analytic  in  the  cut  plane  (cut  from  0  to  oo  along  the  real 
axis).  A(z)  does  not  have  this  property,  but  we  may  write 


u>0 

A-(«)  X-(u) 


(4.81) 

where  X  (s)  is  analytic  and  non-vanishing  in  the  cut  plane.  An  appropriate 
X  (z)  is 


(4.82) 


Note  that 


Jam  XU)  -/ 

I 


(4.83) 
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In  terms  of  the  X  function,  (4.80)  becomes 


E(o)ia.[cf^u)  ^ ^  «  X V)aV)-x'(u)a*(u)  u>o 


(4.  84) 


Thus  J  (z)  defined  by 


85) 


is  an  entire  function.  But  since  we  require  that  C((z)  vanishes  as  z  — >co, 
J  (z)  must  also  vanish  at  oo.  Hence 


(4.  8M 


and 


.  J _ I  r”x-m  [cF(u)-9.^“jh4^,. 

(X(i)  x(i)  •nci  A'(“)  u-a 

with  the  desired  expansion  coefficient  A  (  v)  given  by 


(4.87) 


(4.  88) 


The  normalization  constant  C  may  be  obtained  by  letting  z  — roo  in  (4.87) 


IU> 

HKe 


(4.  89) 


identity 


The  integrals  in  (4.  89)  may  be  simplified  through  the  use  of  the 
18 


XU) -I  =  -L  f . 

ZTCl  J  K  -Z 

o 

_  j_  r‘°X'(u)[A*(u)-A'(u)1  I  _ 

'  A^Wfw-z)  ~"w‘-jAru64-a) 


(4.90) 
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which  follows  from  the  fact  that  X  (z)  -  1  is  analytic  in  the  cut  plane, 
vanishing  as  |z|--^  oo.  For  our  purposes,  we  need  only 


X(o)- 


-«E  f  uF'(u)dM 

x'w .  - 

<*>  /,  A'N 


Then,  using  the  definition  of  g.(u)  (4.  38)  and  the  fact  that 


aaF(vi)  = 


(4.91) 


(4.  92) 


(4.  93) 


we  may  obtain  the  normalization  constant 

which  could  be  further  simplified  slightly  since  we  know  (cf.  eq.  (  3.67]) 


x(o)  =  (I  - 


(4.  95) 


7.  Conclusion  -  General  Considerations 

We  have  carried  out  the  eigenfunction  expansion  and  its  application  to 
boundary  value  problems  in  plasma  oscillations.  It  is  seen  that  the  equations 
treated  contain  a  continuous  spectrum  of  allowed  values  with  associated  singular 
normal  modes.  In  other,  more  fundamental  respects  -  completeness  and  ortho¬ 
gonality  to  the  accompanying  normal  modes  of  the  adjoint  problem  -  the  eigen¬ 
function  expansion  for  plasma  problems  does  not  differ  from  the  eigenfunction 
expansions  that  arise  in  all  branches  of  physics.  Thus  the  main  effort  in  this 


1 


dissertation  has  been  the  attempt  to  place  boundary  value  problems  in  plasma 
oscillations  on  an  equal  footing  with  other,  perhaps  more  familiar,  boundary 
value  problems. 

The  existence  of  a  continuous  spectrum  with  singular  normal  modes  is 

not  peculiar  to  the  plasma  problem.  Actually  it  is  a  characteristic  outgrowth  of 
D 

the  -  or  "streaming"  operator  that  appears  in  the  equations  governing  a 

t 

24 

wide  variety  of  transport  phenomena.  Along  this  line.  Case  and  Dyson 

have  presented  the  continuum  modes  for  the  Euler  equations  of  hydrodynamics 

25 

(linearized  about  a  linear  velocity  profile).  Also,  Cercignani  has  obtained 
the  solutions  to  certain  boundary  value  problems  in  gas  dynamics  through  an 


expansion  in  the  (singular)  normal  modes  of  the  linearized  Boltzmann  equation 
for  a  gas  after  approximating  the  scattering  kernel  by  a  kernel  of  finite  rank. 
Finally,  Case^  has  provided  a  comprehensive  singular-normal-mode  treatment 

li 

of  one-velocity  neutron  transport  theory,  and  his  results  have  been  applied  and 


generalized. 


27-30 


In  Case's  treatment,  the  time -dependent  homogeneous  neutron  transport 


equation  with  plane  symmetry  and  isotropic  scattering 


(4.  96) 


under  the  ansatz 


(4.  97) 


takes  on  the  familiar  form 


(4.  98) 
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uisUncet  are  measured  in  unit  of  the  absorption  mean  free  path,  time 
in  units  of  /v  and  c  is  the  average  number  of  neutrons  emerging  from  a 
scattering  event. 

The  solutions  to  (4.  98)  are 

a)  for  any  k,  a  continuum  -  1  ^  ^  1 

-  I  + 


(4.  99) 


b)  for  |k(  one  discrete  mode 


'f'  zk  «;-/t 

v/ith  OC  the  root  of  A{oc^  )  =  0  or. 


(4.  100) 


i  -  Un(i(^) 


(4.  101) 


The  orthogonality  relatione  are 


-I  *  ' 

^  z 

Jk* 


(4.  102) 

The  norn^al  mode  expansion  may  be  used  to  obtain  the  general  time 
dependent  Green's  function  -  the  neutron  distribution  due  to  a  source 


.108. 


q  =  S  (t  -  t  J  5  (x  -  X  )  S  (  M  -M  ).  For  our  purposes,  it  suffices  to  state 

O  O  /  *  0  I 

a  less  general  result:  the  total  neutron  density  ^(x,  t)  s  produced 

by  an  isotropic  plane  source  q  =  $  (t  -  t^)  $  given  simply  by 

^  •iff  <  ^•f.k  (4.  103) 

(here  X  means  adding  the  discrete  term  to  an  integral  over  the  continuum). 

« 

There  is  an  interesting  connection  between  this  neutron  transport 
problem  and  the  transport  problem  of  the  gaseous  discharge.'^  Let  'Y\v,  v^,  z,  t) 
be  the  distribution  function  for  electrons  in  an  electric  field  diffusing  among  gas 
molecules.  The  electric  field  is  in  the  z  direction  and  plane  symmetry  and 
elastic  isotropic  scattering  is  assumed.  Then  >|/*  satisfies 


(4. 104) 


where  and  (T  are  the  absorption  and  scattering  cross-sections.  If, 

•  s 

finally,  we  assume  that  the  cross-sections  are  proportional  to  1/v  and  choose 
units  in  terms  of  the  mean  free  time  (v  (X"^)  and  the  basic  acceleration  eE/m, 


we  get 


-  Y  +  T  Y 


/ rdK  % 


(4. 105) 


where  c  =  CT  /  (T  . 

s  a 


We  seek  a  Green's  function  solution  q  =  S(v)  S  (z)  ^(t)  corresponding 
to  the  injection  of  a  pulse  of  electrons  at  rest.  Since  the  solution  will  contain 
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a  factor  S(v  - by  conaervation  of  energy;  no  information  ia  loat  by 
integrating  (4.  105)  over  z  to  obtain 

^  St  •  r  4  ^  ^ 


(  4. 106) 


where 


107) 


-»  •  (4. 

Now  if  we  use  instead  of  v,  v^,  the  velocity  variables  v  and  yL4  =  v^/v 
the  transport  equation  takes  on  the  form 

I  i  =  ^ 

which  is  identical,  under  the  transcription  v — •.r  to  the  neutron  transport, 
equation  in  spherical  geometry.  Thus,  the  solution  $  corresponding  to  a 
point  source  in  velocity  space  is  obtained  inomediately  by  transcription  from 
the  neutron  transport  point  source  Green's  function.  Thus,  the  results  of 
the  eigenfunction  expansion  in  neutron  transport  may  be  applied  directly  to  the 
gaseous  discharge  problem. 
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